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Abstract. The prototypical problem we study here is the following. Given & 2L x2L square, 
there are approximately exp(4i4'L^/7r) ways to tile it with dominos, i.e. with horizontal or 
vertical 2x1 rectangles, where K « 0.916 is Catalan's constant |12l I22|. An algorithmically 
simple way of sampling uniformly one among so many tilings is to introduce a Markov Chain 
algorithm where, with rate 1, two adjacent horizontal dominos are flipped to vertical dominos, or 
vice- versa. The unique invariant measure is the uniform one and a classical question |23l 1141 113j 
is to estimate the time Tmix it takes to approach equilibrium (i.e. the running time of the 
algorithm). In [14) . fast mixing was proven: T^ix = 0(L'^) for some finite C. Here, we go much 
beyond and we show that c\L^ < Tmix < C2L^^''- Our result applies to rather general domain 
shapes (not just the 2L x 2L square), provided that the typical height function associated to the 
tiling is macroscopically planar in the large L limit, under the uniform measure (this is the case 
for instance for the Temper ley- type boundary conditions considered in U). Also, our method 
extends to some other types of tilings of the plane, for instance the tilings associated to dimer 
coverings of the hexagon or square-hexagon lattices. 



1. Introduction 

Uniformly random perfect matchings (or dimer coverings) of a bipartite, infinite, planar peri- 
odic graph G (e.g. 7? or the hexagonal lattice Ti) play a crucial role in statistical mechanics and 
combinatorics, and a vast literature exists on the subject (cf. for instance the classical papers 
\12\ [22] and the much more recent [8]). On one hand, thanks to the bijection between perfect 
matchings and discrete height functions (see Section [2. they provide natural and exactly solv- 
able models of random (2 -|- l)-dimensional interfaces (which can be thought of as simplified 
models for the interface separating two coexisting thermodynamic phases L20j). On the other 
hand, thanks to their conformal invariance and Gaussian Free Field- like fluctuations properties 
in the scaling limit [911101(8]. they belong, like the Ising model at T = Tc, to the family of critical 
two-dimensional systems. 

In contrast, the study of stochastic dynamics of perfect matchings is a much less developed, 
though equally interesting, topic. Typically, one takes a large but finite portion G' of the graph 
G and defines a simple Glauber-type Markov chain such that each update modifies locally the 
matching within G' and the unique equilibrium measure is the uniform measure over perfect 
matchings of G' . From the point of view of theoretical computer science |14| [131 [23l I17j . the 
interesting question is to understand how quickly, as a function of the size of G' , the Markov 
chain approaches equilibrium (i.e. how quickly it samples reliably a uniformly chosen random 
matching of G'). From the point of view of statistical physics, thanks to the above mentioned 
bijection between matchings and height functions, this Markov chain can be seen as a dynamics 
for a (2 -|- l)-dimensional interface and it is of interest to understand how the geometry of the 
interface evolves in time. For instance, when G = % then the evolution of the height function 
exactly coincides with the zero-temperature heat-bath dynamics for an interface, separating "— " 
from "-|-" spins, of the three-dimensional nearest-neighbor Ising model [1]. 
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Until recently, mathematical results on this dynamics issues (in particular [13]) were limited 
to proving that, when G = 7? oi G = T-L, the equilibration time Tmix of the Markov chain is 
at most polynomial in the size of G' . These results are based on simple but effective coupling 
arguments that unfortunately have little chance of providing the sharp behavior of Tmix- (A 
notable exception is the work [23], where sharp bounds on Tmix for G = % were given, but for 
a very particular, spatially non-local, dynamics). 

In [5], instead, in the case where G = % it was proven that, under a certain condition 
on the shape of the finite region G' ("planar boundary condition" assumption, see below), 
Tmix behaves like L^, up to logarithmic corrections, if L is the diameter of G' . The major 
conceptual improvement with respect to [HI US] is that [S] crucially uses the intuition that 
the height function should evolve, on diffusive time-space scale, according to a deterministic, 
mean-curvature type evolution [20]. By the way, such intuition suggests the precise scaling 
Tmix ~ const X L? log L. 



As we briefly explain in Section 1.2, going beyond the case of the hexagonal lattice ?^ is a 
real mathematical challenge, since certain "miracles" that occur there do not survive on more 
general graphs. In this work, we prove that, when G = 1? (and in a few other cases, see below), 
again under the "planar boundary condition" assumption, Tmix ~ L^- We present this result, 
still informally but with a bit more of detail, in the next section. 

1.1. Informal presentation of the main result. A domino tiling of the plane is a covering of 
with 2x1 non-overlapping vertical or horizontal rectangles (dominos), with vertices sitting 
at points of {1?)* = 1? + (1/2,1/2). Domino tilings are in one-to-one correspondence with 
perfect matchings oi G = 7?, i.e. subsets of edges (called dimers) of 1? such that each vertex is 
contained in exactly one dimer (to see the correspondence, just draw a segment of unit length 
inside each domino, parallel to its longer side, with endpoints on Z^). Similarly, tilings of a 
finite portion P of the plane correspond to perfect matchings of a finite subset G' of Z^. From 
now on, we will abandon the tiling language and adopt the perfect matching one. Typically, if 
the set G' does admit matchings (an obvious necessary condition is that its cardinality is even) 
and its area is large, the number Z{G') of matchings grows like the exponential of a constant 
times its area. This is for instance the case when G' = {1, . . . , 2L} x {1, . . . , 2L}, in which case 
logZ(G') ~ AKL'^/t:, where K ^ 0.916 is Catalan's constant [T21I22]. 

A natural way to uniformly sample one among so many matchings is to run a Markov chain 
where, with unit rate, two vertical dimers belonging to the same square face of n G' are 
flipped to vertical, or vice-versa. The unique stationary (and reversible) measure being the 
uniform measure over all matchings of C, a classical question in theoretical computer science 
\14:\ [T3l [23] is to evaluate how quickly the Markov chain reaches equilibrium, as a function of 
the diameter (call it L) of G' . This is measured for instance via the so-called mixing time Tmix 
(see Section [5] for a precise definition). 

The only known mathematical result we are aware of on this problem is in ^14j: Tmix = O(-L^) 
for some C < 00 (no lower bound was given). In this case, in the theoretical computer science 
language, the Markov chain is said to be "rapidly mixing" (slow mixing would correspond to a 
Tmix growing faster than polynomially in L). 

Our result in the present work is much more refined and we show that the correct constant is 
C = 2. The price we have to pay is a non-trivial restriction ("planar boundary height" condition) 
on the shape of the region G' , that we briefly introduce now. The lattice being a bipartite 



graph, it is possible to associate in a canonical way (see Section 2.1) a discrete height function 
(defined on faces of G') to each perfect matching of G' . The height along the boundary dG' of 
G' is instead independent of the matching and depends only on the shape of G' . We will say 
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that the boundary height of G' is planar if the graph of the height function, restricted to dG' , 
is within distance of order 1 from some plane of M^. In this case, for L large the height function 
of a typical matching of G' (sampled from the uniform measure) is macroscopically planar not 



only along dG' but also in the interior of G' (see Theorem 3.7). 

The planar boundary height hypothesis is not so exotic and it is verified for instance when 
G' is the 2L x 2L square as above. More general domain shapes that verify this hypothesis 
are introduced in [9] ( "Temperley boundary conditions" ) and in that case the height function is 
proven to converge to the Gaussian Free Field |Hl HH] . 

Our main result can be informally stated as follows (see Sections 2.1 and [5] for a precise 
statement of the hypothesis and of the result): 

Theorem 1.1. // the diameter of G' is L and the boundary height is planar then, as L goes to 
infinity, 

ciL^ ^ Tmix ^ C2(e)L2+^ for every e > 0. (1.1) 

Based on the "mean curvature motion" heuristics mentioned above, we conjecture the true 
behavior to be Tmix ~ const x logL for reasonably regular domains G' . 

We will see that our result holds also when I? is replaced by some other natural infinite 
bipartite periodic graphs, like the hexagon and square-hexagon lattices of Fig. [!} For the 
hexagonal lattice the upper bound in ( |1.1[ ) was known [5], while the best known lower bound 
was Tmix > -^^/(cslogL) (based on [3]). We are not aware of previous results for the square- 
hexagon lattice. 



As we explain in Section 9.3, however, there are good reasons why we cannot consider general 
bipartite periodic planar graphs (for instance, why our method necessarily fails for the square- 
octagon graph of Fig. [l]). This is related to the existence for such graphs of so-called "gaseous 
phases" in their phase diagram [8j. In a gaseous phase, the height function looks qualitatively 
like a (2 -|- l)-dimensional low temperature Solid-on-Solid interface (the interface is rigid, height 
fluctuations have bounded variance and are their spatial correlations decay exponentially. In 
the scaling limit, the interface does not behave like the Gaussian Free Field in this case). 

Remark 1.2. The main reason why in Theorem \l.l\ we require the boundary conditions to be 
almost-planar is that in this case the height fluctuations at equilibrium (i.e. under the uniform 
measure) are well- controlled, see Theorems \3.(}\ and\3.'7\ In the case of general boundary condi- 



tions, only partial results are known (e.g. lll| i and these are not sufficient to implement our 
scheme. We emphasize that instead the Tmix = 0{L'-") result of |14j does not require boundary 
conditions to be almost-planar. 



1.2. Sketch of the proof and novelty. Here we briefly sketch how the proof of Theorem 1.1 
works, and we point out the main novelties, especially with respect to flM 15]. 

The idea of [H |S] is to break the proof of the upper bound Tmix < cL?'^'^ into two steps: 

(i) first prove that Tmix ^ cL?'^'' when the height function is constrained for all times between 
a "floor" and a "ceiling" that are at small mutual distance, say L"^/^"; 

(ii) then, via an iterative procedure that mimics the mean curvature motion that should emerge 
in the diffusive limit, deduce Tmix < cL?'^'^ for the unconstrained dynamics. 

While this general scheme is robust and will be employed also here, step (i) is very much 
model-dependent. In particular, in [H [5] for the hexagonal graph % its implementation was 
based on the crucial observation (by D. Wilson [23]) that, for a certain non-local version of the 
dynamics introduced in [13], one can write explicitly an eigenfunction of the generator, and the 
evolution of the height function is controlled by the discrete heat equation. As we explain in 
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Remark 8.8, this miracle fails for graphs other than T-L and it has to be replaced by a more 
robust argument. 

The "more robust argument" starts from the observation that, under the non-local dynamics, 
the mutual volume Vt between two evolving height functions, that can be seen as an integer- 
valued random walk, decreases on average with time t. This was already realized in [14j, but 
without any further input this only implies a polynomial upper bound ©(L^"^"^) on the mixing 
time (the maximal volume between two configurations in a region of diameter L, constrained 
between floor and ceiling at mutual distance U/^^ , is of order L^"*""^/^", and a random walk 
Vt with negative drift will hit the value with high probability in a time at most of order 
(^2+e/i0)2 < ^4+E_ ^j^gj^ ^j^g volume is zero, the two configurations have coalesced and a simple 
coupling argument allows to conclude that Tmix < L^^*^). The new input we provide (see Section 



9.1 ) is that the variance of an increment (V^+i — Vt) can be lower hounded essentially by Vt itself: 
then, a martingale argument allows to conclude that the coalescence time is of the (optimal) 
order L^"'"'^ and not L^"*"*^. 

To prove the bounds on drift and variance of Vt, we introduce a mapping between perfect 
matchings and configurations of what we call a "bead model". This mapping turns out to be 
very convenient in that it makes the proofs visually clear. In particular, the definition of the 
non-local dynamics is much more natural in this language than in the "non-intersecting-path" 
language of |il4j . 

A last comment concerns the mixing time lower bound in Theorem which is better (by a 
factor logL) than the lower bound Tmix > CL^/logL found in |5]. First of all, the proof of [5] 
would not extend for instance to G = Z^, again because it is based on Wilson's eigenfunction 
argument that fails there. Moreover, even for G = % where Wilson's argument works, removing 



the log L in the denominator involves a genuinely new idea, see Section 10 : one needs to prove 
that the drift of the volume Vj, which as we mentioned is negative, is not smaller than the size 
of the boundary of G', times some negative constant. 

1.3. Organization of the paper. All the definitions and results we need about perfect match- 
ings, height functions and translation-invariant infinite measures of given slope are in Sections 
[2] to |4] (here results are given for general periodic bipartite planar graphs and not just for the 
square, hexagon and square- hexagon graphs). The dynamics is precisely defined in Section [5] 
and its monotonicity properties are discussed in Section [6| In Section [7| we map height functions 
into the configurations of what we call a "bead model" . In Section [8] we rewrite the dynamics in 
terms of beads and we introduce two auxiliary, spatially non-local, dynamics, that are essential 
in proving the mixing time estimates of Theorem |1.1[ the mixing time upper bound is proven 
in Section [9] and the lower bound in Section [Tol 



2. DiMER COVERINGS, HEIGHT FUNGTIONS AND UNIFORM MEASURES 

We follow the notations of |8]. Let G = {E,V) be an infinite, Z^-periodic, bipartite planar 
graph. "Bipartite" means that its vertices can be colored black or white in such a way that 
white vertices have only black neighbors and vice-versa. Z^-Periodicity means that G can be 
embedded in the plane in such a way that Z^ acts as a color-preserving isomorphism. The dual 
graph of G, whose vertices are the faces / of G, is denoted G* . 

We let Gi (the fundamental domain) denote G/1?, which is a finite and periodic bipartite 
graph, embedded on the two-dimensional torus. See Fig. [l] for some classical examples (the 
square, hexagon, square-octagon and square-hexagon lattice) together with their fundamental 
domains. 
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Figure 1. Some examples of Z^-periodic bipartite graphs (square, hexagon, 
square-octagon, square-hexagon) and their Newton polygons. The fundamental 
domain is indicated with thicker lines in the graph while the action of is 
represented by two arrows. Near each vertex of the Newton polygon is indicated 
the associated matching of the fundamental domain. 

Note that there is a certain degree of arbitrariness in the embedding of G in the plane and as 
a consequence a certain arbitrariness in the choice of the fundamental domain. For instance, in 
Fig. [ijthe fundamental domain of contains two sites, but with a different embedding it could 
contain the four sites sites around a face (in this case the two axes of Z^ would be horizontal and 
vertical). In general, it is convenient to work with the smallest possible fundamental domain, as 
in Fig. [TJ 

A perfect matching (or simply "matching" in the following) of G is a subset of edges, M G E, 
such that each vertex of G is contained in one and exactly one edge in M. It is known that 
G admits a matching (which is implicitly assumed from now on) if and only if Gi does, and 
Fig. [T] shows that the fundamental domains of the four graphs we mentioned do admit several 
matchings. We denote the set of matchings of G. 

Assumption 1. To avoid trivialities, we will assume that for every edge e of G there exists 
M £ Q such that e £ M and M' G i7 such that e ^ M' (one can easily construct pathological 
examples where this fails, hut the edges in question can he simply removed and the matching 
prohlem is unchanged). 

We will often refer to paths on the dual graph G*: 

Definition 2.1. A path 7 on G* is a possihly infinite sequence (. . . , /o, /i, . . . ) of faces of 
G, such that fi is a neighbor of /j+i . An infinite path 7 is called periodic if there exists a finite 
path 70 and v £ I? such that 7 is the concatenation 0/ {Tn^7o}nez, with Ty the translation by v. 

2.1. Height function and uniform measure. A flux is a function on the oriented edges of 
G, which is antisymmetric under the change of orientation of the edges. To each M G is 
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associated a flux um- edges contained in M carry unit flux, oriented from the white to the black 
vertex. Edges not contained in M carry zero flux. Note that the divergence of ujm is 1 at white 
vertices and —1 at black vertices. 

Fix now a reference matching Mq (typically, a Z'^-periodic matching, but the following def- 
inition would work for any flux of divergence +1/ — 1 at white/black vertices). Mq allows to 
associate to M a height function Hm on G* , as follows. Fix some face /o G G* ("the origin") 
and set /iAf(/o) to some value, say 0. For every / ^ /o, let 7 be a path on G* starting at /o 
and ending at /. Then, /ia/ (/) — hMifo) is the total flux of lom — ^Mq (say from right to left) 
across 7. Note that Km does not depend on the choice of the path (because ojm — comq has zero 
divergence) and that the height difference between two matchings M, M' is independent of the 
choice of the reference matching Mq. 

In the following, the set of perfect matchings will denote equivalently the set of all admissible 
height functions (it is understood that Mq and /q are fixed). For lightness of notation, we will 
often write h instead of Km- 

Definition 2.2. Let U be a simply connected open subset set of [—1/2, 1/2]^, L > and Ul = 
L U . We let G' be the finite subset of G obtained by keeping all the vertices and edges belonging 
to faces which are entirely contained in Ul- 

Given m £ il. (called the "boundary condition", with corresponding height function h^) and 
G' = {E',V') as in Definition [2^ let 

n^,G' = {M £ n : M\g\g' = m\G\G'} (2.1) 

be the finite collection of matchings that coincide with m outside of G' . Equivalently, we can 
identify r2m,G' as the set of height functions that coincide with km except on the faces of G' . 
We will implicitly assume (without loss of generality) that the reference face /o is not one of the 
faces of G' . Clearly, i^rn,G' is non-empty (it includes at least m) and we will let Tim^G' denote 
the uniform measure over ^rn,G' ■ 

3. Pure phases 

In this section we review known results about measures on the infinite graph G whose typical 
height functions are close to a plane. First we will identify the set of "natural" measures of fixed 
average slope, then give a classification into three "phases" with very different correlation prop- 
erties, and finally give their "microscopic" behavior, i.e. the probabilities of events depending 
on a finite subset of edges. Most results come from [8] . 

3.1. Ergodic Gibbs measure of fixed slope. Fix a reference matching Mq assumed to be 
Z^-periodic. A measure on is said to have slope (s,t) if its expected height function is 
linear of slope (s,t): for all faces /, if /' denotes the translate of / by (x,y) G 1?, then 
fi[h{f') — h{f)] = sx + ty. fi is said to be a Gibbs measure if its conditional distributions on 
finite subgraphs are uniform, fi{ ■ |M G ^m,G') = ^m,G'(') {DLR property). It is ergodic if it is 
not a linear combination of other Gibbs measures. Ergodic Gibbs measures of fixed slope can 
be thought of as the natural uniform measures on matchings of G conditioned on their average 
slopes. The following theorem due to Sheffield jl8j classifies all of them : 

Theorem 3.1. There exists a closed, convex polygon N in such that, for all {s,t) in its 
interior N, there exists a unique ergodic Gibbs measure fis.t of slope {s,t). The vertices of 
N are determined by the slopes of some 'Z'^ -periodic matchings of G (i.e. matchings of the 
fundamental domain) and thus are integer points. For {s,t) £ dN := N\ N, there exists an 
ergodic Gibbs measure but it may not be unique. 
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is called the Newton polygon, see Fig. [T] 

3.2. Phase classification. As proved in [8j, ergodic Gibbs measures come in three possible 
phases: solid, liquid and gas, depending on the position of (s, t) in A^. 

• Solid phases correspond to slopes in dN. For any side I of dN, there exists at least 



one infinite periodic path 7 on G* (cf. Definition 2.1) such that the configuration of the 
edges crossed by 7 (or by any of its translates) is deterministic and is the same for all 
measures with slope (s, t) G I. The path 7 is said to be frozen. 

The asymptotic direction of 7 is determined as follows. All the planes with slope in 
i and containing the origin of intersect in a straight line. The direction of this line, 
when projected on the (x, y) plane, is the direction of 7. 

At a vertex of the Newton polygon, which is the intersection of two sides of dN , there 
are two families of frozen paths with different directions, which form so to speak a grid 
on G* . The components of the complement of the frozen paths are finite sets of faces. 
Heights are clearly independent in two distinct components, and the fluctuations of the 
height difference between two faces /i , /2 are bounded deterministically and uniformly 
in the distance between them. 

Liquid phases correspond to generic points of iV- In these phases, heights fluctuations 
behave like a Gaussian free field in the plane. In particular the variance of the height 
difference between /i and /2 grows like l/vr times the logarithm of the distance, while 
edge correlations decay slowly (as the inverse of the square of the distance). Liquid 
phases are discussed in finer detail in the next section. 

Gaseous phases have exponentially decreasing edge correlations; the height difference 
fluctuations are not deterministically bounded, but their variance is bounded, uniformly 
with the distance of the faces. Gaseous phases may (but do not necessarily) occur when 
the slope (s, t) is an integer point in N . The condition fo r th e occurrence of a gaseous 
phase at an integer slope (s,t) GiV is discussed in Section 



3.4 



In the example of Fig. [T| only the square-octagon graph has a gaseous phase which has slope 
(0,0). 

3.3. Edge probabilities. When (s,t) GiV (i.e. for liquid and gaseous phases) there is an 
explicit expression of edge probabilities under ij,s,t- 

Theorem 3.2. [8] Fix {s,t) SiV- There exists an infinite periodic matrix Ks^t = {-f^s,t(w, b)}w,b 
with b (resp. w) ranging on black (resp. white) vertices of G and an infinite periodic ma- 
trix K~l = {K~l{h,w)}^ Ta satisfying Ks^tK~l = Id such that, for any finite subset {ei = 
(wi, bi), . . . ,ei = (w;, bz)} of edges of G, the Hs,t-pi"obability of seeing all of them occupied is: 

^ls,t{el e M, . . . , Q E M) = ^J] A',,t(wj, b,)^ det(i^-i(bfe, w,)), ^ . ^ ^ ^. 

Kg^t is called a Kasteleyn matrix. It is a weighted and signed version of the adjacency matrix, 
so in particular i^<j^t(w,b) can be non-zero only if (b,w) is an edge of G. The signs (which are 
independent of the slope {s,t)) are chosen so that their product around any face / is — 1 if / 
has mod 4 sides and +1 if it has 2 mod 4 sides. We will not need to specify the explicit 
choice of signs, see 0. Periodicity means that Ks^t{^ + {x, y),h + {x, y)) = -fCs^t(w, b) for every 
{x,y) E Z^, and similarly for K~l . 

Given Kg^t and two complex numbers w, z, we define a finite matrix K.g^tiz, w) from white to 
black vertices of the fundamental domain Gi, as follows. Consider Gi as a weighted periodic 
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bipartite graph on the torus, where the weight of an edge is the one induced by Ks^t, and note 
that it can contain multiple edges between two vertices, even if the infinite graph G does not 
(see e.g. Fig. [T]). Consider a path '~^x (resp. 7y) winding once horizontally (resp. vertically) 
along the torus and multiply by z (resp. 1/z) the weight of each edge crossed by 7^ with the 
black vertex on the left (resp. on the right) and similarly by w,l/w the edges crossed by 7^. 
Then, 'Ks,t{z,w) is the adjacency matrix of Gi, with these modified weights. With the usual 
graph theory convention, this means that the (w, b) element of 'K.s,t{z,w) (with w (resp. b) a 
white (resp. black) vertex of Gi) is the sum of the weights of the edges joining w to b. Let 
Qs,t{z,w) (a matrix from black to white vertices of Gi) be the adjugate matrix of K.s^t{z,w) 
so that [Qs^t'i^s,t]{z,w) = P{z,w)ld where P{z,w) = det{'Ks,tiz,w)). The (b, w) element of 
Qs,t{z,w) is denoted Q^'^ {z,w). 

We can now give a formula for the inverse infinite Kasteleyn matrix K~l: 

Theorem 3.3. [8j Let h and w be a black and a white vertex in Gi. The following holds for 
(x,y) GZ2.- 

T^-i^u , ^ ^^ 1 /" '5s,f(^>w^) X ydwdz 
(2z7r)^ Jr2 Ps,t{z,w) w z 

where = {z, w G (D^, \z\ = \w\ = 1} is the unit complex torus. 

Here, to avoid confusion, it can be useful to emphasize that [Qs,t/ Ps,t]{z,w) is the inverse 
of the finite matrix 'Ks,t{z,w), while K~l is an inverse of the infinite matrix Ks,t, with no 
dependence on z,w (i.e. with the original edge weights). 

3.4. Asymptotics of K~l and Gaussian fluctuations in the "liquid phase". It is shown 
in [8] that, for an integer slope {s,t) £N CiT?, the Laurent polynomial Pg^t has either no zeros 
on the unit torus (in which case fis,t corresponds to a gaseous phase) or has a unique zero of 
order two (which corresponds to a liquid phase). For any non-integer slopes in instead, 
Pg^t has exactly two conjugate simple zeros. In this case Theorem 3.4 gives the asymptotics of 
K~l{h,w) when the two vertices b,w are far apart. We emphasize that in our applications (i.e. 
in the proof of Theorem 3.6), we will have to consider only cases where the Newton polygon has 
no integer points in its interior. 

Theorem 3.4. [8j Fix {s,t) a non-integer slope in N , so that P^.t has two simple zeros {zo,wq) 
and {zo,wo) on T^. Let a = -g^Ps^t{zo,wo) and /3 = -§^Ps,t{zo,WQ) and define (j){x,y) = xazQ — 
y^WQ. Then the map (f> : IR^ 1— >• (D is invertible, the matrix Qs,t{zo,wo) is of rank 1 and can be 
written as Ug^ty^t where the column vector Ug.t (resp. Vs^t) is indexed by the black (resp. white) 
vertices of Gi . Moreover, we have 

«.-(.,w.,„,).-.(«M^).o(^) 

where 0((x^ + y'^)~^) has to be understood as ^^i^y^^^ with h bounded on 1? and 9(2;) denotes 
the imaginary part of z. 

Remark 3.5. The invertibility of (p is a consequence of the fact that azQ is not colinear with 
Pwq. This is not proved explicitly in [8] but the argument is simple: Both the torus and 
P^tii^})' set of zeros of Pg^t, (^fc two-dimensional manifolds in (D^ which contain {zq^wq). 
The tangent space ofT'^ at {zq,wo) is given by 

pi = {a{izo, 0) + 6(0, iwo) + {zo,wo), (a, b) G M^} 
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and the tangent space of Pg^{{0}) is 

P2 = {C{P,-a) + {zo,wo),C G(D}. 

Since {zo,wo) is a simple zero of Pg^t seen as a function on T^, one necessarily has pi r\p2 
{{zojWq)} and it is easy to check that this fails if ozq = \(3wq with A € M. 



The asymptotic expression (3.1) is the main tool for the fohowing result that will be proved 
in Appendix [A} 

Theorem 3.6. Fix a non-integer slope {s,t) GN and g G G* . Under the Gibhs measure iJLs,t, 
the moments of the variable 

h{f) - h{g) - ifisAHf)) - f^sAHg))) h{f)-h{g)-{ii,AKf))-l^s,m9))) .3 

^Var,^AKf)-h{g)) ^\ogm)-m\ 

tend as \(j){f) — — t- 00 to those of a standard Gaussian J\f{0, 1). 

Here and later, when we write (/>(/) we mean (j){x, y) if / is the (x, y) translate of a face in the 
fundamental domain Gi. The fact that the variance of h{f) — h{g) behaves like (I/tt) log !</>(/) — 
4>{g)\ is proved in [8]. 

For the hexagonal lattice and under the assumption that / and /o are along the same column 
of hexagons, convergence of the moments is proven in [7J. As we explain at the beginning of 
Appendix [AJ the general case is qualitatively more difficult and requires non-trivial work; our 
proof uses ideas from [llj Sec. 7]. 

3.5. Almost-planar boundary conditions. A central role will be played by "almost planar" 
boundary conditions. 

We say that /i G is an almost-planar height function with slope {s,t) £N if there exists C 
such that, for every f £ G* , 

\hU)-h^s,t{Km<C. (3.3) 

We will sketch briefly in Section [4] a proof that almost-planar boundary conditions actually 
exist for every GiV (even with C = 1). 
Theorem 3.6 implies the following: 

Theorem 3.7. Fix a non-integer slope {s,t) GN and let m be an almost-planar boundary 



condition with slope {s,t). Let the finite graph G' be as in Definition 2.2. One has for every 
e > and n > 0; 

vr„,G'(3/ e G* : |/i(/) - /i.,t(/i(/))| > L') = 0{L-^). (3.4) 

Proof of Theorem \3. 7| given Theorem \3.6\ For the hexagonal lattice this is given in detail in [4J 
(see Proposition 4 there). For general graphs the proof is almost identical and we recall just the 
basic principle. 

By monotonicity (see Section [6]), the event Ef = {h{f) — ids,t{h{f)) > L'^} is more likely 
if we change the boundary condition for a higher one, i.e. if we replace m with m' such that 
hm'if) ^ hfnif) for every face /' ^ (G')* adjacent to some face in (G')* (this set of faces is 
denoted here dG' , and {G')* is the collection of faces of G'). Assume without loss of generality 
that the reference face /o where heights functions are fixed to zero belongs to dG' . Then choose a 
random boundary condition m' from the measure Hs,t, and this time fix its height at the reference 
face as /im'(/o) = ^m(/o) + L^f^ = Thanks to Theorem 3.6, one has hm' > hm on dG' , 

except with probability 0{L^^) for any given n. Finally, with such random boundary condition, 
by the DLR property the probability of Ef is nothing but ^s,t{hif) ~ l^s,t{h{f)) > /2), which 
is also 0{L~'^), again thanks to Theorem 3.6 □ 
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4. Perfect matchings, capacities and maximal configurations 

4.1. Linear characterization of height functions. The set of height functions corresponding 
to a perfect matching of a finite subset of G can be characterized by linear inequalities as fohows. 

Consider as in Definition |2.2| a finite subgraph G' of G and a boundary condition m G fi. In 
this subsection we will use m (even if it is not necessarily periodic) as reference matching for the 
definition of height functions. For any two neighboring faces /, /' with common edge e oriented 
positively (i.e. such that going from / to /' one crosses e leaving the white vertex on the right), 
let the oriented capacities d{f, f) and d{f', /) be defined as follows: 

'Oife^G' fOife^G' 
d{f, /) = < if e G G' and e G m d{f', /) = < 1 if e G G' and e G m . (4.1) 

1 if e G G' and e ^ m; if e G G' and e ^ m 

Now for any pair of faces /, /' (not necessarily neighbors) let D(f, /') be the minimum over all 
paths / = /i, ...,/„ = /' in G* of the sum of the d(/j, fi+i) (the minimum is well defined, the 
capacities being nonnegative) . 

Proposition 4.1. An integer-valued function h on G* is the height function (with reference 
matching m) of a matching in ^ra,G' "if O'l^d only if 

Dif, f) ^ hif) - h{f) for every f, f e G* . (4.2) 
Proof. The proof is in the spirit of [6, Theorem 1]. The "only if" part is trivial since, going back 



to Section 2.1, it is immediate to see that the maximal possible height difference h{f') — h(f) 



between neighboring faces /, /', for any matching in fi^.G') does not exceed d{f,f'). As for 



the "if" part, remark first of all that, thanks to (4.2), for every neighboring faces /, /' one has 
h{f) — h{f') G {—1,0, 1}. Let us "mark" all edges e G G' between faces /, /' (with e oriented 
positively from / to /') such that /i(/') — h{f) = d{f, /'), together with edges e £ G \ G' such 
that e G m. Let M be the union of all marked edges and let us prove it is a matching (note 
that, automatically, M = m outside of G'). For any white (resp. black) vertex v, let e„ be 



the unique edge incident to v which belongs to m. From (4.2) and considering paths that turn 
anticlockwise (resp. clockwise) around v, it is easy to see that: 

• either all the faces / sharing vertex v have the same value of h(f) and e„ is the single 
marked edge around v; 

• or there exists a single marked edge e'^ ^ e^, incident to v, such that h{f') = h{f) — 1, 
with /, /' neighboring faces sharing e^, such that v is on the left (resp. right) when going 
from / to /'. 

M is thus a matching with by constructionM = m outside G'. In conclusion, M G ^m,G' 

of course its height function is just h. □ 

4.2. Maximal and minimal configurations. The characterization of height functions pro- 
vided by Proposition |4.1| shows the existence of a unique maximal (resp. minimal) height function 
hmax (resp. hmin) in ^m,G'- "Maximal" means that for any other height function h in ^rn,G' 
satisfying /i(/o) = (recall from Section [s] that the height is fixed to zero at some face /o outside 
of G') one has h{f) < hmax{f) for every f £ G* . Indeed, define hmaxif) ■= D{fo, f) on G*. This 



satisfies (4.2) (since D{-,-) satisfies the triangular inequality) and maximality is a consequence 
of the fact that d{f, f) is the maximal possible height difference between neighboring faces. 
Similarly, one has hminif) = —D{f,fo). Observe that the height functions hmax, hmin (with 
respect to the reference configuration m) vanish outside G' as they should (this is because the 



set of faces of G not belonging to G is connected, recall Definition 2.2). 



DYNAMICS OF PERFECT MATCHINGS 



11 



4.3. Free paths and possible rotations. 

Definition 4.2. Fix a matching M G fi. We say that an oriented path 7 in G* is a free path 
(relative to M) if all edges crossed by 7 are free (i.e. not occupied) and have the same orientation 
(i.e. either all of them have their white vertex on the right of j or all of them on the left). If 
white vertices are on the right (resp. left) then 7 is called a positive (resp. negative) free path. 

See Fig. [2} 




Figure 2. The relation between height function and dimer configuration in the 
visually appealing case of the honeycomb graph. It is easy to see the right-hand 
drawing as a stepped surface in 3 dimensions. It should be clear from the right- 
hand drawing that the positive free path (thick line) moves away at constant 
speed from the (1, 1, 1) plane; at its endpoint (marked by an arrow) the free path 
cannot be possibly continued, and a cube can be added there (i.e. a rotation can 
be performed in the left-hand drawing). 

A first observation is that free paths cannot form loops: 

Proposition 4.3. Let ^ be a free path relative to some M G il, and assume that 7 forms a 
sim,ple loop. Then, for every M' S $7 the edges crossed by 7 are free. 

Together with Assumption [T| this excludes loops. 

Proof of Proposition Let h be the height function of M', with reference matching M. Let 
/ be a face along 7. By symmetry, suppose 7 is a positive path. Since all edges are traversed 
with the positive orientation, we have 

h{f) = h{f) + Hedges crossed by 7 and occupied in M'}| 
— Hedges crossed by 7 and occupied in M}|. 

Since 7 crosses no occupied edge of M by assumption, it crosses no occupied edge of M' either. 

□ 

A second observation is that, since only the reference matching (however it is chosen) makes 
a contribution to the height difference along a free path 7, the height function is non- increasing 
(resp. non-decreasing) if 7 is a positive (resp. negative) free path. An important consequence, 
that we will need in Section [8] to upper bound the equilibration time of the dynamics, is the 
following: 
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Proposition 4.4. Fix {s,t) £N- Let M £ Q be such that the corresponding height function 
stays between two planes of slope {s,t) and mutual distance H. All free paths relative to M have 
length at most CH, where the constant C depends only on {s,t). 



Proof of Proposition 4-4 ■ Since the graph G is periodic, there exists only a finite number k of 
types of faces that are not obtained by integer translation of each other. Let 7 be a positive free 
path (if it is a negative free path, the argument is similar) relative to some matching M. We 
claim that, 

if one walks n steps along 7, the function / 1— h{f) — fis,t{h{f)) decreases by at least (4.3) 

— [n/fcje for some e = es,t > 0. 

Then, the proposition follows (with the constant C being inversely proportional to e/k) because 
the function ^s,t{h{-)) on G* is essentially planar with slope {s,t). 

To prove ( |4.3| ), observe first that the matching M gives no contribution to the variation of 
h along 7 (all crossed edges are free) so that the variation of /i — ^s,t{h) is simply minus the 
/Xs,t-average number of crossed edges which are covered by dimers. Fix some face / G 7 and walk 
along 7 until a face /' which is a translate of / is reached (the number of steps is at most k). 
The /Ug^f-average of crossed dimers between / and /' is non-negative and we will actually prove 
that it is strictly positive and independent of the type of face /, which implies the claim. Indeed, 
let 7 be the infinite periodic path on G* obtained by repeating periodically the finite portion of 
7 which joins / to /'. If the average of crossed edges is zero, then clearly the slope of the height 
under the measure ^s,t along the asymptotic direction of 7 is extremal, which contradicts the 
assumption that (s, t) is in the interior of the Newton polygon A^. Uniformity w.r.t. the type of 
the face / is just a consequence of the fact that the number of different face types is finite. □ 

For any face / of G, there exist exactly two ways to perfectly match its vertices among 
themselves. Label "+" one of the two matchings, and "— " the other (according to some arbitrary 
rule). If M is a matching of G such that the vertices of / are matched only among themselves, 
we call "rotation around /" the transformation which consists in leaving M unchanged outside 
of /, and in flipping from "— " to "+" (or vice-versa) the matching of the edges of /. If some 
vertices of / are matched to vertices not belonging to /, then the rotation is not possible. 

Free paths yield a way to flnd a face where an elementary rotation is possible. Given M G fi, 
we pick an arbitrary face /i and we construct a growing sequence {7n}n>i = {(/ii ■ ■ ■ ■, fn)}n>i, 
of positive free paths, with 71 = (/i) (an analogous construction gives a growing sequence of 
negative free paths). Given 7„, consider all faces / which are neighbors of /„ and such that going 
from fn to / one crosses a free edge with white vertex on the right. Choose fn+i (according 
to some arbitrary rule) among such faces. If there are no such faces available, we say that the 
procedure stops at step n. In this case, it means that every second edge around fn is occupied 
by a dimer, and this is exactly the condition so that a rotation at fn is possible. Altogether, we 
have proven: 

Proposition 4.5. Fix {s,t) GiV- Let M £ Q be such that the corresponding height function 
stays between two planes of slope (s,t) and mutual distance H. Within distance Gg^tH from any 
face f there exists a face f~^ (resp. f~ ) where a rotation is possible; such a rotation increases 
(resp. decreases) the height at (resp. f^) by 1 and h{f^) ^ h{f) ^ h{f~). 

4.4. Almost planar height functions. Here we prove that almost-planar height functions 



satisfying (3.3) with C = 1 do exist (under the assumption that {s,t) GiV is a non-integer 



slope). Indeed from Theorem 3.6 and Borel-Cantelli we get that, for every fixed 5 > 0, almost 
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all configurations from fis,t satisfy 

m) - i^sAHm <B+ mf) - (4.4) 

for some random B, where /o is the face where the heights are fixed to zero. Take one of these 
configurations. Let An the set of faces at graph-distance at most n from /o and suppose that 
h{f) — (j,s^t{h{f)) < —1 (the argument is similar if the difference is > 1) for some / G An- 



The same argument that led to Proposition 4.4 shows that, if 5 is chosen small enough (say 



much smaller than the constant Sg^t in (4.3)), any positive free path 7 starting from / is of 
length 0{5n/ £s^t) ^ n for n large enough. Therefore, the last face /' of 7 is in A2n and (by the 
properties of positive free paths) one has h{f') — ^s,t{h{f')) < —1. By Proposition 4.5 , a rotation 



is possible at /' and it increases h{f') — ij,s,t{h{f')) by 1. The configuration thus obtained clearly 



still verifies (4.4) with the same B and the quantity 



A= ^ \h{f) - ^ls,t{h{f))\l\h(f)-^,,^,(h{f))\>^ 

decreased by 1. Since A is finite, the procedure can be repeated a finite number of times until 
there is no point left in An with \h{f) — fis,tih{f))\ > 1. One concludes easily using the fact that 
n can be taken arbitrarily large. □ 

5. Dynamics and mixing time 
The dynamics we consider lives on the set ^}rn,G' of matchings on a finite subset G' C G (as in 



Definition 2.2 ) with boundary condition m G 0. Every face / of G' has a mean-one, independent 
Poisson clock. When the clock at / rings, if the rotation around / is allowed, flip a fair coin: if 
"head" then choose the "+" matching of the edges of /, if "tail" then choose the "— " matching. 
In other words, perform the rotation around / with probability 1/2. 
Call /i*^ the law of the dynamics at time t, started from M. 

Proposition 5.1. For t — )■ 00, fj,^^ converges to the uniform measure vr^^G'- 

Proof. It is obvious that TTm,G' is invariant and reversible, so one should only check that the 
dynamics connects all the configurations in O^.g'- This is done by using the free paths of 
Section HI 

Let M G ^m,G' and let M^ax £ ^m,G" be the matching corresponding to the maximal height 



function /imax introduced in Section 4.2 The height function h oi M with reference matching 
Mmax is clearly non-positive and vanishes outside G' . Pick a face / such as h{f) < —1 and 
consider a positive free path 7 growing from / (as in the proof of Proposition |4.5| ) . Along 7 the 
height function h cannot grow and, G' being finite, 7 has to stop after a finite number of steps. 
The last face /' of 7 clearly is inside G' (since the height is zero outside) and we have already 
discussed that a rotation is possible at /' and it increases h{f') by 1. By recursion, M can be 
transformed into Mmax by a finite sequence of elementary rotations inside G' . Arbitrariness of 
M allows to conclude. □ 

As usual [13j, an informative way to quantify the speed of approach to equilibrium is via the 
mixing time, defined as 

Tmix = T;nix(m, G') = inf{t > : max Wfif - vr^.c'll < l/(2e)} (5.1) 

where — is the total variation distance of measures /x, and the choice of the value l/(2e) 
is conventional (any other value smaller than 1/2 would do). With this choice, one has [13j 

max 1 1 /if -irmG'\\<e~ WTmi^i . (5.2) 
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We will study the mixing time when the boundary 
conditions are "almost planar". The following is the 
main result of this work: 

Theorem 5.2. Fix {s,t) £N, m an almost-planar 
boundary condition of slope {s,t) and let G' be as in 
Definition \2.<^ If G is either the square, hexagon or 
square-hexagon lattice (cf. Fig. [7]j then there exists 
some ci > and, for every e > 0, C2 < oo such that 



ClL < Tmix < C2L 



2+e 



(5.3) 



For the hexagonal lattice the upper bound was 
known |5j, while the best known lower bound was 
Tmix ^ -^^^/(cs log L) We emphasize that removing 
the logarithm in the denominator involves a genuinely 
new idea. For the square lattice, it was only known 
that Tmix = 0{L^) for some C < 00 jH] but no lower 
bound was given. We are not aware of previous results 
for the square-hexagon lattice. 



Figure 3. An example of the class 
of graphs where our results could 
be extended, see Remark |5.4[ The 
shaded region is a thread (cf. Sec- 
tion [7]). Layers of squares and 
hexagons can be of arbitrary verti- 
cal thickness and the periodicity in 
the "vertical" direction can be arbi- 
trarily large. 



Remark 5.3. The proof of the lower bound in (5.3) 



actually shows the following: if the dynamics is started 
from the maximal configuration, which has an excess 
volume cL/^ with respect to the typical (almost flat) 
equilibrium configuration, it takes a time ciL^ before 
the excess volume becomes smaller than say (c/2)L^ 
(which is still very large w.r.t. typical volume fluctua- 
tions). In this sense, the equilibration time lower bound 
is optimal. 



Remark 5.4. Our result could be extended to a class 
of graphs obtained by alternating periodically layers of 
squares and hexagons (see Fig. On the other hand, we will explain in Section 9.3 why our 
method does not ( and should not!) work for general periodic bipartite graphs G, in particular 
not for graphs like the square- octagon lattice which possesses a "gaseous phase". 

6. MONOTONICITY 

It is natural to introduce the following partial order on i7 : M ^ Af if and only if /ijv/(/) ^ hM'if) 
for every f £ G* . As usual, the reference face /o is assumed to be fixed once and for all. Note 
that the partial order does not depend on the reference matching used to define the height. We 
say that an event A C is increasing if M ^ M' and M' G A implies M £ A. We define in the 
usual way stochastic domination: fi ^ fi' if /u(j4) ^ fJ-'iA) for every increasing event A. 

Proposition 6.1. The dynamics defined in Section^is monotone, that is iif^ ^ /x^^' for every 
t and every M > M' . 

Proof. Couple the dynamics started from M and M' by using the same clocks and the same 
coin tosses. Partial order is preserved along time. Indeed, it suffices to observe that if /im(/) = 
hM'if) ^ rotation at / that increases the height by 1 is possible for M' , then necessarily the 
configuration of the edges of / in M is the same as in M' , otherwise at some face /' neighboring 
/ one would have /ia/(/') < hM'if')- ^ 
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Remark 6.2. As in O Sec. 2.2], one can realize all the evolutions for all possible initial 

conditions Mq on the same probability space, with the property that if Mq < Mq then almost 

surely M/^« < Aff ° for every t > 0. This construction is called global monotone coupling. 
Proposition 6.3. If A is an increasing event, then TTm,G'{'\A) >z iTm,G'- 



Proof. Remark that in the proof of Proposition 5.1 we showed that the maximal configuration 



can be reached from any other by a chain of rotations that increase the height, so /imax £ ^ and 
A is connected. Consider the original dynamics started from /imax and the reflected dynamics 
(again started from /imax) where each update that would leave A is canceled. It is clear that 
they converge to tt^^g' and T^m,G'{'\-^) respectively and that, when coupled by using the same 
clocks and coin tosses, the second always dominates the first. □ 

Monotonicity allows to apply "censoring inequalities" of Peres and Winkler [15] which, roughly 
speaking, say the following: if the dynamics is started from the maximal or minimal configu- 
ration, deleting some updates along the evolution in a pre-assigned way (i.e. independently of 
the actual realization of the dynamics) increases the variation distance from equilibrium. The 
precise statement we need (cf. Corollary |6.5| below) is a bit more general than what is proven 
in [15] but the proof is almost identical, so we will just point out where some modification is 
needed. 

Consider a probability measure vr on 17 and P^'"'^ , w G V a set of transition kernels that satisfy 
reversibility {'K{a)P^'"\a rj) = ■K{rf)P^'"\r} — )■ a)) and monotonicity. We define a dynamics on 
O by assigning a Poisson clock of rate c„ to each f G V and applying p(^) when v rings. The 
dynamics of Section [s] corresponds to V = {G')* , = 1 and P'^^^ the kernel that corresponds to 
a rotation around / with probability 1/2 (if allowed). 

Theorem 6.4. Let z/q be a probability measure on Q, such that ^ is increasing. Consider vt 
the law at time t of the dynamics started from vq. Then, for every t >0, ^ is increasing and, 
if {fit}t>o is a family of probability measures such that vt ^ for all t, one has 

\\vt - 7r|| ^ ll/it - 7r||. 



Corollary 6.5. Let uq be as in Theorem 6.4- Suppose that for all u G V, for all v such that 
^ is increasing, we have P^'"^v < v. Let /xj be the law at time t of the dynamics started from 
vq, where the rates Cy of the Poisson clocks are replaced by deterministic time- dependent rates 
Cv{s), such that < Cv{s) < c^ for every < s < t. Then, 

for every t > 0, ft ^ l^t o-nd \\ut — 7r|| ^ \\^t — 

Remark 6.6. The hypothesis o/di^o/dvr increasing is immediate if the dynamics is started from 
the maximal configuration h^g^^, since in that case is concentrated on hmax- 



As in [15], the proof of Theorem 6.4 follows directly from the following two lemmas. 



Lemma 6.7. With the above definitions, for any probability measure n, if ^ is increasing then 



IS increasing. 



This replaces Lemma 2.1 of [15j . which uses explicitly the fact that the dynamics is of "heat- 
bath" type. 
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Proof. 



(^) 



1 
1 

Tr(a) 



1 

vrfcr 



^ 7r(s) 



5:glvr(.)p(^)(.^.)=E(^)[^(X)] 



(6.1) 
(6.2) 



where X is the state after one action of P^^\ starting from a. The third equahty uses the 
reversibihty and the monotonicity of P^^-* shows that the last expression is increasing in u. □ 

Lemma 6.8. [151 Lemma 2.4] If ii,v are two probability measures on 17 such that ^ is increasing 
and y ^ jJL, then \\v — tt\\ ^ — 7r||. 

Proof of Corollary \6.5\ Decompose the Poisson point process (PPP) of density on as the 
union of two independent PPPs, X and Y, of non-constant densities Cv{s) and — Cv{s) > 0. 
The dynamics nt is obtained by erasing the updates from the process Y. From Lemma [6. 7| we get 
that ^ is increasing. Censoring an update P^^^ at time t conserves the stochastic domination 
because, by induction, ut = P^^^u^- ^ u^- ^ /i^- = fXf 



□ 



7. From dimers to a "beads'' 



From this point onward, we will assume that the graph G is either the square, hexagon or 
square-hexagon graph (see Fig. [T]) since we will use some of their geometric properties. 

The set Q of matchings of G can be mapped into the configurations of what we call a bead 
model. Such correspondence is valid for more general graphs than the square, hexagon and 
square-hexagon, provided that the graph in question possesses a certain "fibration" with fibers 



(or threads) satisfying the properties described below. See Remark 5.4 for a more general class 
of graphs where this construction would work. 




Figure 4. The paths called threads and used in the definition of the bead model 
for the square, hexagon and square hexagon graph. The arrow shows the orien- 
tation of the threads. 

As is apparent from Fig. [4j for the three types of graph we are considering, there exists a 
family of directed periodic paths {7i}iGZ (called threads) on G* such that 
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(i) labeling the faces along thread 7^ as {/*}jgz, face fjj neighbors only and some faces 

of7i±i; 

(ii) going from the face /j to /j+i, one crosses an edge of G (call it e* ) which is positively 
orientated. Such edges are called transversal edges and dimers on transversal edges are 
called beads; 

(iii) threads 7^ are obtained one from the other by a suitable 7? translation and Uj7i = G*. 



Proposition 7.1. Consider a finite subgraph G' C G as in Definition 2.2 (recall that G\G' is 
connected) and a boundary condition m £ 0,. A matching in ^m,G' ^-^ uniquely determined by 
the position of its beads. Furthermore the number of beads in G' on each thread ji is the same 
for every M £ il.m,G'- 

Proof. From the definition of height function and using property (ii) above of threads, which says 
that transversal edges are all crossed with the same orientation, the height at some / G {G')* 
belonging to thread 7^ is determined by the number of beads on 7j between the boundary of 
G' and /. Hence the position of beads uniquely determines the height function and thus the 
matching. The total number of beads on 7jn(G')* is determined by the height difference between 
two faces /i, /2 of 7, \ (G')*, such that the portion of 7, between fi and /2 includes 7^ n {G')* . 
This height difference is clearly independent of the particular chosen matching in fim.G' (because 
fi , /2 are outside G' and G \ G' is connected) . □ 

In order to have a complete picture, we have to determine the condition for a set of beads' 
positions to correspond to a matching of G, that is the kind of constraints beads impose on each 
other. We first need some notations (see Fig. [s]). 

Definition 7.2. Let as above e*- denote the j*^ transversal edge of the thread 7^ and let b*- (resp. 
Wj) denote its black (resp. white) vertex. Let F* be the set of vertices of G which belong both to 
7j and to 7j+i, and order vertices in F* following the same direction as for the faces along the 
threads. Note thatV^ contains the white vertices of transversal edges ofji and the black vertices 
of transversal edges of ^i-i^i. Given transversal edges e^^^e^^^ on 7^,74+1, we write e^j ^P'^ (resp 
e* >- e^j1~^) ifw^j is below (resp. above) b^-^^. 

Proposition 7.3. A set of bead positions corresponds to a matching in 0, if and only if, for any 
two consecutive beads on the same thread (i.e. beads on transversal edges e\,e\.,a < b with no 
bead between them along the same thread ji), there is a unique bead in thread ■ji-i and a unique 
bead in thread ^i+i such that their positions e*~^, e^^^ satisfy e\ -< e*~^ -< e\ and e\ -< e!"^^ -< e^. 

Proof. We advise the reader to keep an eye on Fig. [5] while reading this proof. 

Proof of the "only if" part. Without loss of generality we can consider threads 7^ and 7^+1. 
Note that following F* between w^ and w^ there is necessarily exactly one more black vertex than 
white vertex (because G is bipartite) . Since by assumption there are no beads on 7, between ejj 
and el, all the white vertices have to be matched within F*. This leaves exactly a single black 
vertex which has to be matched along a transversal edge in 7i+i. The corresponding dimer is 
the unique bead such that its position e^^^ satisfies -< e^"^^ -< e^. 

Proof of the "if" part. Suppose that the bead positions are given and that they satisfy the 
properties above. This automatically fixes which transversal edges are occupied and which are 
free. To see that the rest of the matching is also (uniquely) determined, proceed as follows. 
With the same notations as above, consider the vertices of F* between w^ and b^"^^. These 
vertices are all matched with each other (because by construction e^^^ is the "highest" bead in 
7j+i such that e^^^ -< el) and they form a path with an equal number of alternating white and 
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black vertices, so there is a unique way of matching them. The same goes for vertices between 
w^ and hi+\ □ 

Remark 7.4. Fix G' C G and a boundary condition m £ 
Under the measure 7rm,G' ; conditionally on the positions of the 
beads in 7j±i, the beads of are independent and each has a 
uniform distribution in a certain finite set of adjacent transver- 
sal edges ( two transversal edges being adjacent if they are of the 
form e*-, e*_,_]^; remark that on the square lattice they can actu- 
ally share a vertex, while on the hexagonal lattice they cannot). 
Indeed, given a bead on the transversal edge e\, it is possible to 
move it up to e\_^i (resp. down to e\_i) via a rotation of the 
face fl^i (resp. provided that fl^^ (resp. /*J is within G' 
and that the new position does not violate the ordering properties 
of Proposition [TT^ Uniformity of the distribution is trivial from 
uniformity of the unconditional measure TTm,G' ■ Note also that 
moving a bead up (resp. down) implies changing by —1 (resp. 
+1) the height of the face just above (below) it along the thread. 

8. Dynamics in terms of beads 

The Glauber dynamics defined in Section [5] has a simple in- 
terpretation in terms of bead model. As observed in Remark |7.4[ , a rotation is equivalent to 
moving a bead to an adjacent transversal edge in the same thread (in particular, rotations are 
possible only at faces adjacent to a bead, and the configuration of beads outside G' is frozen). 
We can then redefine the dynamics as follows. Each bead in G' has a mean-one independent 
Poisson clock; when it rings, with probabilities 1/2, 1/2 move the bead either up or down to the 
adjacent transversal edge if allowed by the boundary conditions and by the ordering properties. 

8.1. Fast dynamics. As mentioned in the introduction, we will not work directly with the 
original Glauber dynamics but rather with an auxiliary one. We will actually need two auxiliary 
dynamics: one, that we call synchronous fast dynamics, will be useful to upper bound the mixing 
time, while the asynchronous fast dynamics will provide a lower bound. 

Definition 8.1. 

(i) We define the synchronous fast dynamics as follows. We have two independent mean-1 
Poisson clocks. When the first (resp. second) one rings we resample all the bead posi- 
tions on even-labelled (resp. odd-labelled) threads following the equilibrium measure 7rm,G' 
conditioned on the state of beads on odd (resp. even) threads. 

(a) The asynchronous fast dynamics is defined instead by giving each bead in G' an independent 
mean-1 Poisson clock. When a clock rings, the position of the corresponding bead is resam- 
pled from iTm^G' conditioned on the position of all other beads. Note that, by construction, 
beads outside G' are frozen. 

Thanks to monotonicity of the original Glauber dynamics, one sees easily that both synchro- 
nous and asynchronous fast dynamics are monotone. 

Remark 8.2. Recall Remark \ 7.4\ the positions accessible to a single bead, given the beads of 
neighboring threads and the boundary conditions, form a segment of the transversal edges of its 
thread and these segments are mutually non-intersecting. Thus, under both synchronous and 
asynchronous dynamics each bead is resampled with the uniform law on a finite segment. 




Figure 5. An example of 
the situation in the proof of 



proposition 7.3 Beads on the 



edges e\ and are indicated 
by thick lines. 
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8.2. Comparisons of mixing times, and constrained dynamics. As announced, the asyn- 
chronous dynamics will provide a mixing time lower bound for the original one. 



Proposition 8.3. Fix G' <Z G as in Definition 2.2 and a boundary condition m, not necessarily 



almost planar. Let Tmix he the mixing time for the original dynamics and let T" he the first 
time such that the law of the asynchronous dynamics, started from the maximal configuration, 
is within variation distance l/(2e) from CQuilihrium, Then, T^[^ 



Proof. Recall the definition (5.1) of mixing time: to get a lower bound, we can just look at the 
evolution from the maximal configuration hmax, which we can therefore assume to be the initial 
configuration of both original and asynchronous fast dynamics. From the description of the 



asynchronous fast dynamics in Section 8.1, we see that we can couple the two dynamics using 



the same Poisson clocks for each bead. To prove that the asynchronous dynamics approaches 
equilibrium faster than the original one it is enough to show that if ^ is increasing then, writing 
P^^^ and P^'^^ for the kernels corresponding to an update of the bead b according to the original 
and asynchronous fast dynamics respectively, one has 

pWzy ^ pWz.. (8.1) 

Indeed, together with monotonicity this guarantees that the height function is stochastically 
lower under the asynchronous dynamics than under the original one and then Theorem |6.4| can 
be applied (with fit the law of the original dynamics and vt that of the asynchronous one). 
By conditioning on all the beads except 6, we can assume that is a measure on some interval 



{0, . . . ,k} and that vr is the uniform measure on the same interval (cf. Remark 7.4) in which 



case it i s tr ivial to check (8.1). Indeed, pC'V = vr for every u, and dP^^V/dvr is increasing (cf. 



Lemma 6.7 and recall the assumption di^/dyr increasing), which implies vr ^ P^^^v. □ 



Constrained dynamics. Next, we bound Tjam from above using the synchronous dynamics: this 
works well only if the dynamics is constrained between two configurations whose height functions 
are not too different. Given two matchings M_ < M+ in ^rn,G' (M_ will be called "the floor" 
and M+ "the ceiling" ) the constrained dynamics is defined in the subset VL^ q, C J^m,G' such 
that M_ < Af < M+: it is obtained from the original dynamics, erasing all updates which would 
exit ^^^^Qi- It is elementary to check that monotonicity still holds, and the equilibrium measure 

is of course tt^"^, := 7rm,G'("l^- ^ " ^ ^'-^+)- The distance between floor and ceiling is defined 
as H = maxjgG-. (/ij\,f^(/) — hM_{f))- To avoid a proliferation of notations, W6 still Cctll T^mix 

the 

mixing time of the dynamics constrained between M_ and M_|_, and fif^ its law at time t, with 
initial condition M. 

To estimate Tmix within logarithmic multiplicative errors, we can restrict ourselves to the 
evolution started from the extreme configurations (see for instance [H Eq. (6.5)]): 

Lemma 8.4. Consider the dynamics constrained between M_ and M+. For any t > and any 



l^f -^^G'IK2i?|(G')*|max 



I M+ _ M± II II M_ _ M± 
\l^t ■^m,G' II' II /^i '^m,G' 



'm,G' 

with |(G')*| the number of faces of G' . 

In the next result, T^- denotes the mixing time of the synchronous dynamics constrained 



between M_ and M_|_ (just take the Definition 8.1 of the synchronous dynamics and replace 
■^m,G' by iT^%, there): 
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Proposition 8.5. Fix G' C G as in Definition 2.2 and boundary condition m. Suppose m is 
almost planar with slope {s,t) £N and consider floor /ceiling M± G Qrn,G' > c-t distance H from 
each other. We have 

T^i. ^ {C,,tH'log^H\{Gr\))T^i.- 

Proof. A complete proof for the hexagonal graph is given in [4] Section 6.2] and since for the 
square or square- hexagon graph not much changes, we will be somewhat sketchy. 

For simplicity we will write vr for tt^%i. One first proves that a single update of the syn- 
chronous fast dynamics can be realized by letting the original dynamics evolve for a time 
0{H'^ (log H){log\{G')*\)) while censoring some updates. Indeed, consider the dynamics ob- 
tained by starting from M and setting to the rate of updates of beads on, say, odd threads in 
the original dynamics. It is clear that this auxiliary evolution converges to the uniform measure 
on configurations of beads on even thread conditioned by the beads on odd threads, which is 



exactly the measure after one "even update" of the synchronous fast dynamics, and Remark 8.2 



allows us to easily compute its mixing time. Beads on odd threads are frozen and those on even 
threads are completely independent so we have to compute the mixing time for a set of indepen- 
dent one- dimensional simple random walks on domains of the type {0, . . . , ki}. By Proposition 
4.4 the ki are bounded by Cg^tH (because, if a bead can be moved n steps up, necessarily there 



is a length-n free path along the thread) so the mixing time for each walk is 0{H^ log H) and 
for 0(|(G')1) such walks it becomes 0{H^ {log H){log\{G')*\)). 

It is then clear that the law of the synchronous fast dynamics at time t, call it v^^, coincides 
(except for a negligible total variation error term), with that of the original dynamics at time 
n{t)A after censoring suitable updates, where 

A = Cs,tHHlogH){log\{Gr\) < Cs,tH^log\H\{Gr\) 

and n(t) is a Poisson random variable of average 2t (this is the number of updates within 
time t for the synchronous dynamics). Since we will take t large, we can replace n(i) with 
its average (we skip details) . If Corollary |6.5| is applicable (see below) , we obtain then that 



M± 



vr < Wv, 



M± 
t 



tt\\. Then, using Lemma 8.4 and (5.2) 



II M 
SUp||^2AAT= 
M 



max 



\l^2AATt 



-7r|| < 2H\{G' 

X 

< 2H\{G' 

which is smaller than l /(2e ) for some A of order log(i^|(G') 



I M_ 



max 



^AT" 



^AT" 



< 2H\{G' 



(8.2) 
(8.3) 



To see that Corollary 



6.5 



is applicable, we have to check that d/i/dvr increasing implies P^^^ fj, ^ 
fi with P^'') the kernel of the update of bead b under the original dynamics, where beads move 
by ±1 along their respective thread. Conditioning on all other beads, we can assume that n is 
a probability on an interval {0, . . . ,k} and that vr is the uniform measure on the same interval. 
Then, summation by parts shows that 

1 

x=l 



which is negative if / is increasing {fi is also increasing). 



□ 



8.3. Volume drift. In this section we study the time evolution of the volume between two 
configurations under the (a) synchronous fast dynamics. This will be the key to evaluate their 



mixing time and thus, thanks to Propositions 8.3 and 8.5 the mixing time of the original 
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dynamics. Note that in Proposition 8.6 we do not require the boundary condition to be almost- 
planar. 

Proposition 8.6. Let M^,M^ £ nm,G' be such that ^ and let Ml,i = 1,2 denote the 
evolution starting from and following the fast dynamics (synchronous or asynchronous: we 
use the same notation). Letting Tt denote the filtration induced by {M],i = l,2,s < t} and 
Vt = J2 feG*i^M^(f) ~ ^A/i(/)]; then Vt is a supermartingale: 

E[Vt\Tt'] ^ Vt> for every t > t' . 
The same holds if the fast dynamics is constrained between a floor M_ and a ceiling M+. 

Note that, since the volume is expressed as a sum of height differences, it does not matter 
whether M} evolves independently of Mf or not. 

Proof. First of all note that the expected drift 

\un-mVe)-V,] (8.4) 

e\0 e 

is the same for the synchronous and for the asynchronous fast dynamics (as a function of the 
initial conditions M^,M^): this is thanks to the fact that the volume is a sum of height dif- 
ferences over faces, that expectation is linear and that beads in the same thread are updated 
independently in a step of the synchronous dynamics. (This does not imply that the process Vt 
itself or even its average is the same for the two dynamics) . In the following, we will therefore 



assume that we deal with the synchronous dynamics and prove that (8.4) is not positive. 



From the proof of Proposition 5.1, we see that there exists a sequence of configurations 
M(Q), . . . , M(fc) such that M(o) = M'^, M^^) = and M(j) is obtained by M(j_i-) via a rotation 
that decreases the height at some face /. Writing the volume difference between and 
as a telescopic sum of volume differences between M(j-) and M(j_i) and using linearity of the 
expectation, we see that to prove (8.4) we can restrict to the case where and differ only 



by a rotation on a single face /. We will actually prove that the expected change of V from a 
single update is except when / is suitably close to the boundary of C, in which case it can be 



negative (see Remark 8.7 for a more precise discussion). 

The proof will be given only for the square- hexagon graph since it contains all the difficulties, 
but the same method works equally well for the hexagon or square graph. 

Case 1: / is a hexagon. Suppose that the face /, whose rotation brings to M^, is a 
hexagon on thread 7j (cf. Fig. [g]). Then, there is a certain bead b which in the two configurations 
is on two different adjacent transversal edges of 7^ on the boundary of / (such edges are called 
ei and 62 in the picture). Now consider a step of the synchronous dynamics. If threads with 
the same parity as 7^ are updated then clearly the evolved configurations can be coupled in 
order to coincide, and the volume decreases by —1. We need to show that when threads of the 
opposite parity are updated the average volume increase is at most +1. Clearly, only threads 
7j_l-i contribute to the average change of volume (because there is no discrepancy between 
and on threads 7i±2) which therefore "screen away" ^i±m ?^ > 3 from the discrepancy at /) 
and by symmetry we will just show that 7j_i contributes at most 1/2. 

On 7j_i there are a certain number of beads: call b"*" the lowest bead which is above b (with 



the ordering convention of Definition 7.2) and b~ the highest bead below b. Note that one 
or both of them could be absent in M^, (for instance 7j_i could contain no beads at all): 
however, suitably changing the dimer configuration outside of G' (which has no effect on the 
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Figure 6. Volume drift when / is a hexagon, (a): a schematic representation 
of threads 7^, 7j„i, 7j_2 with the beads b, 6^, b'. (b): a more detailed view of the 
three threads in question, with the hexagonal face / and the transversal edges 
ei, . . . , 67 mentioned in the proof. 



dynamics) we can always assume that such beads do exist (possibly outside G'). Thanks to 
Proposition 7.3, there exists then a unique bead b' in 7i_2, which is lower than 6"'" and higher 
than b~ , see Fig. [6 a). 

A look at Fig. mb) and Proposition 7.3 suffices to convince that only the following two 
mutually exclusive cases can occur (Fig. [6|^a) corresponds to the first one): 

• if b' is at or above transversal edge 63, then 6+ is at or above edge 64. Then, the 
distribution of 6"*" given all the other beads does not depend on whether b is at ei or 
62, and therefore b~^ (and a fortiori beads on ji-i above b'^) gives no contribution to 
the average volume change. As for b~ , instead, its set of possible positions (which forms 
an interval of adjacent transversal edges of 7i-i, recall Remark 7.4), includes exactly 
one edge more (called 67 in the picture) when b is at ei w.r.t when b is at 62- Since 
the position of b~ is uniform among the available positions and since the average of a 
uniform random variable on {0, . . . ,k} is k/2, the average volume change arising from 
bead b~ is (k + l)/2 — k/2 = 1/2 (recall that moving a bead up by one transversal edge 
results in decreasing the height of a face by —1). It is also clear that beads in ji-i and 
below b~ do not contribute to the volume change, since their set of available positions 
is disjoint from that of b~ (cf. Remark 8.2) and does not change when b is moved from 

ei to 62. 

In this discussion, we ignored so far the fact that rotations for the dynamics are allowed 
only in the subgraph G' . For the bead b~ to reach (after the update) its new available 
position 67 starting from the present one, all the necessary elementary rotations along 
thread 7, have to involve faces contained in {G')*. Otherwise, the edge 67 is actually 
not a possible position for b~ : the set of effectively available positions for b~ is then the 
same when b is in ei or 62, so the volume change associated to 6~ is not 1/2 but 0. In 
conclusion, taking into account the fact that G' is not the whole graph, only decreases 
the average volume change. A similar reasoning shows that the floor/ceiling constraint 
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Figure 7. Volume drift when / is a square. The threads 7j,7i_i,7i_2 with the 
face / and the transversal edges ei,i = 1, . . . , 10. A square face like /, which is 
adjacent to hexagons of the thread to its left, is called a "Type I" square, while 
a square like the shaded one (adjacent to hexagons of the thread to its right) is 
called "Type II". Note that if / were a Type II square say again on 7^, then by 
symmetry ji-i would not contribute to the volume change but the thread ji-^-i 
would. 



can only decrease the average volume change. This observation will be picked up again 
in Remark 18.71 

• if instead b' is at or below transversal edge 65, then b~ is at or below edge sq. The 
argument is then similar to the first case (this time b~ does not feel the effect of the 
discrepancy and 6"*" has one more edge available, again ey, in than in M^). As in the 
first case, constraints from boundary conditions or from floor/ceiling can only decrease 
the volume change. 

Altogether, in both cases the average volume change from 7j_i is either 1/2 or (depending 
on the position of / with respect to the boundary of G' and on the presence of fioor /ceiling) 
and, since the same holds for 7i+i, the overall average volume change is at most zero. 

Case 2: / is a square (cf. Fig. [t]). This time assume that the discrepancy between 
and is at a square face / in 7j, i.e. that a bead b is either on edge ei or 62- Again, when 
threads with the parity of 7^ are updated the average volume decreases by —1. In this case, it is 
also clear that the distribution of beads on ji+i after an update is the same starting from or 
(because the edges ei and 62 meet on the same vertex of 7j+i), so a non-zero contribution 
to the average volume change this time can come only from 7j_i and we will show that this 
contribution is at most +1. With the same conventions as in Case 1 for the beads b^ and b', 
we distinguish this time three cases (to avoid repetitions, we ignore effects due to floor/ceiling 
and to the fact that G' / G: exactly as in Case 1, this can only decrease the average volume 
change) : 
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if y is at or above 63 then is at 64 or higher and does not feel the effect of the 
discrepancy. At the same time, according to whether h is at 62 or ei, the edges 65, eg are 
available or not for h~ . The average volume change induced by h~ is then {k+2)/2—k/2 = 
1. 

similarly, if h' is at or below e^ then does not feel the discrepancy (it is at or below 
eg) and has two extra available positions (again 65,66) when b is at 62- Again this 
gives volume change +1. 

finally, when b' is either at 69 or 610, then both 6+ and b~ feel the discrepancy: indeed, 
when b is at 62 then position 65 is not available for b~ and eg is available for 6+, while 
when b is at 61 then position 65 is available for b~ and eg is not available for 6"'". Again, 
the average volume change is +1. 

□ 



Remark 8.7. It is useful to emphasize that the proof of Proposition 8.6 showed the following. 
Assume for simplicity that there is no floor/ceiling constraint on the dynamics. If and 
differ by a single rotation at f ^ (G')* , then the average volume change after an update is 
between —1 and 0. To decide whether it is zero or non-zero, proceed as follows. The rotation 
around f changes the set of available positions for a certain number (at most two, actually) of 
beads in threads neighboring the thread of f: some positions which were not allowed before the 
rotation of f become allowed and vice-versa. If all the elementary rotations leading such beads 
along their threads to the new available positions are allowed (i.e. if all the faces corresponding 
to such elementary rotations are in the finite domain {G')* ) then the average volume change is 
zero. Otherwise, it is different from zero. 



Remark 8.8. From the proof of Proposition 8.6 one can also understand why the study of the 



dynamics on the square or square-hexagon lattice is qualitatively more challenging than on the 
hexagonal lattice. The basic observation due to D. Wilson in ( although it was not formulated 
in this terms) is the following. Consider the 'fast dynamics" for the hexagonal lattice with initial 
condition given by some matching m, and let IIl^(i) be average at time t of the sum of the heights 
of the faces (in {G')* ) belonging to thread i. Then, for any pair of initial configurations m}, m? , 

1 2 

the function ipt{') ■= HJ" (•) — Hf^ (•) satisfies the discrete heat equation 

dtMi) = l^Mi) + Rt = \[Mi + 1) - '^Mi) + Mi - 1)] + Rt (8.5) 

where the "error term" Rt is due to the boundary of G' and to floor / ceiling constraints (if 
present) and can be ignored for the sake of this discussion. Call Vt = ^i'^t{i) ihe average 
volume difference between the two evolving configurations. Since there are 0{L) threads that 
intersect G' and the lowest eigenvalue of the discrete Laplacian on {1, . . . ,L} is of order L~'^ , 
one deduces immediately that after a time of order x logVo, the average volume Vt is very 
small so that the two configurations have coupled with high probability. 

That the same kind of argument does not work for other lattices can be seen as follows. 
Consider for instance the square-hexagon lattice, and assume that (with the terminology of the 
caption of Fig. and differ only by the rotation at a square face f of Type I on thread 



j. Then, clearly at time zero ipo{j) = 1 and il^oW) = 0, / ^ j. While it is still true that (8.5) 



holds for t = and i = j (the initial drift of'ip{j) is —1 = hA'ip(){j) ), the equation does not hold 



( even at t = 0) for i = j ±1: indeed, the proof of Proposition 8. 6 shows that the initial drift of 
- 1) zs 1 = VV'o(i) = V'o(j) - ^o(j - 1) (instead of 1/2 = \Atpoij - I)) and that of + 1) 
is (instead of 1/2). If on the contrary the face f were a square of Type II, one would find that 
the initial drift of — 1) is and that of + 1) is 1 = —VipoU + !)• 
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We believe that, for initial conditions 'm^,m'^ such that their height differences are "smooth" 



on the macroscopic scale, Equation (8.5) should still (approximately) hold, for L large. Indeed, 
there are as many Type I as Type II squares in each thread: if each of the two types contributes 
approximately equally to the height differences ^poii), from the above reasoning one finds that the 
initial drift ofip{j — 1) is approximately l/2V'i/'o(j) — l/2V'(/'o(i ~ 1) = l/2A?/^o(j — !)■ However, 
trying to pursue this route seems quite hard (one should show that "smoothness" is conserved 
for positive times) and we had to devise an alternative approach instead, based on Proposition 



8.6 and on Theorem \9.1\ of next section. 

9. Mixing time upper bound 



Here we prove the mixing time upper bound of Theorem 5.2 The crucial step (Theorem 



9.1) is to give an almost-optimal estimate when the dynamics is constrained between floor and 
ceiling of small mutual distance H (in the application, we will take H = L'' with e small). Then, 
an argument developed in |5] allows to deduce the mixing time estimate for the unconstrained 



dynamics (see Section 9.2 for a sketch). 



9.1. A martingale argument. The basic step is to prove the following: 



Theorem 9.1. Consider the same setting as in Theorem 5.2, but assume that the height function 
is constrained between ceiling and floor that are almost-planar configurations of slope {s,t) GN, 
of mutual distance H. Then, Tmix = 0{L'^H^log^{HL)). 

The main idea will be to apply to the volume between two configurations the following classical 
bound on the hitting time for a supermartingale: 

Lemma 9.2. Let Xt be a continuous-time supermartingale such that almost surely ^ ^ M 
for every t G and \\m.'mis^Qj¥.[{Xt+s — XtY\Tt] ^ ly > whenever Xt > 0. Suppose 
Xq = i > almost surely, fix < m < i and let T^"^) = inf{t : Xt ^ m} be the hitting time of 
[0,m]. Then we have 

^r„^m^n 2Mi 

E[r(™)] ^ — . 

(Just note that if Z{t) = Xf — 2MXt — min(t, T^''^), then Zt is a negative sub-martingale 
and compute the average of Z{t) for t = T^"*-*). 

Let Vt denote the volume between the maximal and minimal evolutions , Mf under the 



synchronous fast dynamics. Proposition 8^ shows that Vt is a super-martingale. Because of the 
floor and ceiling at distance H, we clearly have ^ ^ \{G')*\H deterministically. To apply 
Lemma 9.2 we only need a lower bound on E[(Vt_|_5 — Vt)^|J-t]. It is important to remark that 
such quantity does depend on how , Mf are coupled, while by linearity it is not necessary 
to specify the coupling to compute the drift E[Vf+5 — Vt|-Ft]. 

Lemma 9.3. There exists a global monotone coupling under which 

liminf lE[{Vt+5 - Vtf\Tt] ^ (9.1) 
where c is a constant depending only on the slope {s,t). 



Proof of Theorem 9.1. Applying Proposition 8.5, it is enough to give the upper bound 

T^i^ ^ cstL^H'^ \og{LH) 
for the mixing time of the synchronous fast dynamics. 
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Let rrii = 2 '^H\{G' )*\ a nd Ti = inf{t : Vt ^ rrii}. Remark that, up to time Tj, Vt satisfies the 
hypothesis of Lemma 9.2 with M = H\{G')*\ and v = crrii/H^. Thus we have 



Finally, since Vt takes integer values, the hitting time of is equal to the hitting time of [0, 1/2] 
which is for iq = [log2(2iJ|(G')*|)l • We have proved 

Er,o ^ c'|(G')*|^^log(i7|(G')*|) = c"L^H''\og{LH). (9.2) 

Therefore, ¥{Ti^ > 2ec" L"^ \og{LH)) < l/(2e), which implies T^;^ < 2ec" L^H'' log{LH): 
indeed, under a global monotone coupling, once maximal and minimal evolutions have coalesced, 
all the evolutions with arbitrary initial conditions have coalesced too. □ 



Proof of Lemma 9.3. Let be the height functions corresponding to the extremal evolutions 
M^. Each face contributes at most H to the volume difference, so there are at least Vt/H 
faces where the height difference hf{f) — h^{f) is at least 1. For each of those, by Proposition 
4.5 there exists a face /~ at distance at most CH where again hf{f~) — h^[f~) > 1 and a 



rotation in that decreases the height is possible. We can thus find at least CVt/H^ distinct 
such faces and each of them is the face directly above a non-frozen bead for (i.e. a bead 
that can be moved upward in via an elementary rotation). Call Bt the set of such beads, 
l^tl > CVt/H\ 

The global monotone coupling mentioned in the claim is defined as follows. We take two mean- 
one independent Poisson clocks: when the first one rings we update the beads in even threads, 
when the second one rings we update beads in odd threads. The beads are updated as follows. 
Suppose for instance that the first clock rings. Then, sample independently for each transversal 
edge e and each bead in each even thread a uniform [0, 1] variable Ub{e) (any continuous law 
would work the same). A bead b in an even thread then chooses the accessible transversal edge 
e (given the positions of beads in odd threads) with the lowest value of Ub{e). It is easy to check 
that this defines a monotone coupling between evolutions with any possible initial condition (we 
emphasize that each evolution uses the same realization of the Ub{e) variables to determine the 
outcome of an update). 

We now turn to the estimate of vt '■= liminf^-^o jIE[(V(4.5 — Vt)^|-F(]. For any bead 6, let vj;^'^ 
denote its contribution to the volume, i.e. the difference of the labels of the transversal edges 
occupied by h in and . Finally let A^~^^ (resp. A^~'>) denote the event that there is an 
update of even parity and no update of odd parity (resp. an update of odd parity and no update 
of even parity) between time t and t + 5 (each has probability 5 + o{5)) and for each bead h let 
A^^'^ be A^^^ according to the parity of h. We have (since the occurrence of two updates has 
probability of order 5"^) 

^n{Vt+5 - Vtf\H = mVt+s - Vtf\J'u A^+^\ + niVt+s - Vtf\J'u A^-^] + o(l) (9.3) 
> Var[l^t+5|Ji,^(+)] + Var[yt+5|J-t,A(-)] +o(l) = ^ Var(y/_;:^| J^, A^^)) + o(l) (9.4) 

b 

(in the last step, we used the fact that conditionally on A^^\ the variables (V^^^ — V^^^) are 
independent for different h and are zero for h of odd parity). For each bead h four cases can 
occur: 

(i) The set of transversal edges accessible to 6 in a single update (given the beads of the other 
parity) is different in and and, at least for one of them, it consists of strictly 
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more than one transversal edge. We let be the set of such beads, 
computatior^ shows that for such bead V'aT:{V^^^g\Tt, A^''^) ^ g. 



An elementary 



(ii) The accessible domain for b is the same in and but its positions in the two 
configurations are different. Let be the set of such beads. Remark that if the event 



occurs, then ^ 



almost surely while V^^^ = V^"' > 1 if it does not. 



(iii) The accessible domain and the initial position of b are the same in 



and Mf- . In this 



case V^_)_^ — "0 conditionally on A^'^^ so these beads give no contribution to the 
volume variation. 

(iv) The accessible domain for b has only a single edge in both and Mf . In this case there 
is no movement possible for b until threads of the opposite parity are updated, so b makes 
again no contribution conditionally on A^''\ 

Remark that the set Bt introduced above is included in BfuB"^, so we have \Bf\+\Bf\^ CVt/H^. 
Indeed, if b ^ Bt than it can be moved in , which excludes case (iv), and / 0, which 
excludes case (iii). 

Suppose that \Bf\ ^ \B^\/{as^tH), with a^.t a slope-dependent constant to be determined 
later; for b G Bf we have Var(VJ^^|J^t, ydC*)) ^ | so the right-hand size of (9.4) gives, taking the 
hmit (5—^0, 



ft ^ 



\B 



7^1 



> a. 



so Vt+s -Vt^ 



IB; 



particular for b G Bf, V^*^^ 



for some positive as^t- 

Suppose on the contrary that \Bf\ ^ \B'f\/{as,tH) and (by symmetry) that at least half the 
beads b G B'f are on even threads. After an even update they each contribute V^^^ — V^''^^ ^ — 1 

2 + Cs^tH\Bf\ ^ — |i?^|/4 (we used the fact that, for every b and in 
^t^^ ^ Cs^t H due to the floor/ceiling at mutual distance H (cf. 
Proposition 4.4) and we chose ag^t = 4cs,j). As a consequence, in the limit 5 — )• (9.3) gives 
i^t ^ bs,tVt^/H^m this case, for some other positive constant bg^t- The conclusion follows from 
mm{Vt/H'^,Vt^/H^) >Vt/H^. ' □ 

Remark 9.4. The power in the lemma is clearly far from optimal. A finer analysis of the 
contribution of Bf and B^ would probably improve the power to . We do not follow this 
route becaus e ul timately the precision of the upper bound is limited by the equilibrium estimate 
of Theorem 3.6 and also because even the bound |-B^| + \Bf\ ^ CV/H^ is certainly far from 
optimal for a typical configuration. 



9.2. A mean curvature motion approach. Given Theorem |9.1| and Theorem 3.7 on the 
equilibrium height fluctuations, the proof of the bound Tmix < c{e)Lp'~^'^ is essentially identical 
to the proof that Tmix = 0(L^(log L)^^) for the dynamics on the hexagonal lattice with almost- 
planar boundary conditions, see |5l Th. 2]. Indeed, Theorem 9.1 plays the role of O Prop. 3] 
while Theorem 3.7 replaces O Th. 1]. Therefore, below we will only recall the main ideas, and 
we skip all details. 



Remark 9.5. The reason why in Theorem 5.2 we get the L'' correction to the mix ing time 
instead of a factor (logL)"*^^ as in [5\ is that the fluctuation estimates of Theorem 3.7 are a bit 
weaker than those of O Th. 1]; since anyway the exponent 12 is certainly non-optimal (we 



^One can check that the worst case is when the intervals of transversal edges accessible to b in and Alf 
are of the form {a, . . . ,a + k — 1} and {a, . . . ,a + k}. In this case, after an update V^''' = with probability 
k/{k + 1) and its average is 1/2 so the variance in question is at least (fc/4)/(fc + 1) > 1/8. 
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conjecture the correct value to be 1, cf. the Introduction), we have not tried to refine Theorem 



3.7 (for instance, one might try to control the exponential moments of the height fluctuations). 

The first step is the following (cf. [5, Prop. 2]): 
Step 1 // the the height function of the initial condition ^ G ^G',m. is within distance U^l"^^ from 
the almost-planar configuration m (i.e. if |^(/) — /im(/)| ^ L"^^^^ for every f G {G')* ), then (for 
any given C ) for all times smaller than the height function stays within distance 2L''^^^ from 
m, except with probability 0{L~'~'). 

This means that, until time , the dynamics is essentially identical to a dynamics with 
floor/ceiling at mutual distance 0{L''^^'^). Together with Theorem 9.1 this implies: 



Step 2 Again if the initial condition ^ is within distance L'^l^^ from m, after time L^"*"^ the law 
of the configuration has small variation distance from equilibrium. 

Therefore, to prove Tmix < c(e)L^+'^ it is sufficient to prove: 

Claim 9.6. At time c'(e)L^"''"^ the evolutions started from maximal/minimal configurations are 
with high probability within distance L''^^^ from m (this is the analog of [5, Prop. 1\). 

Consider for instance the evolution started from the maximal initial configuration hmax £ 
^m,G' and assume for simplicity of exposition that the slope of the quasi-planar boundary 
condition m is (s,t) = (0,0). Let Cu be a spherical cap whose height is u and whose base is a 
disk D of radius pl = L\ogL with the finite graph G' approximately at its center (recall that 
the diameter of G' is of order L <C pi)- Call Ru the radius of curvature of Cu, which satisfies 
{2Ru — u)u = p\ and let ipuif) denote the height of Gu above a face / which is inside D. Then, 
the key to Claim 9.6 is: 



Claim 9.7. With overwhelming probability, the height function (G')* 9 / i— t- /i™"^(/) is below 
the deterministic function {G')* 9 / i— )■ V'm„(/) for all times t G L^], where Un = PL — n and 
the deterministic time sequence tn is defined by 

to = 0, tn = tn-1 + Ru^L'/^, n<M:=pL- 



Indeed, it is easy to verify that Ru„ ~ p\/{2un) and (similarly to [5l Eq. (6.14)]) that 
tM = 0{LP''^''). Therefore, if we show the above claim for n up to M, we deduce that at some 
time 0(-L^+'^) the configuration is within distance um = L''^^'^ from the flat configuration m, 
and Claim [92] follows. 



For n = 0, the statement of Claim 9.7 is true (deterministically, not just with high probability, 
since the maximal height at a face / G G' is of order L <C pi)- Suppose we want to deduce claim 
n+l from claim n, and look for definiteness only at a face / at the center of the disk D. Consider 

1/2 \ 

a disk Dn+i centered at /, of radius Ru^ L with A to be chosen later: by monotonicity, given 
the claim at step n, we can replace the evolution h^^'^^ restricted to -Dn+i, in the time interval 
[fn,i^^], by an evolution where: 

• the configuration outside Dn+i is frozen and equals some height function which, on the 
boundary of -Dn+i, is within distance 0(1) from the function il>u„{-); note that for /' at 
the boundary of Dn+i one has ipunif) ^ Un — {1/2)L'^^'^; 

• the "initial" height function at time tn in Dn+i approximates within 0(1) the function 
V'«„(-)- 

By Step 2, the time such dynamics takes to reach equilibrium is 0{Ru,^L'^^^'^) ^ — tn if 
A < 1/40, so that at time tn+i the configuration is essentially at equilibrium (with the above 
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specified boundary conditions around Dn+i). Next, elementary geometry and Theorem 3.7 
shows that, at equiUbrium, the height function at / is with overwhelming probability lower 
than Un — I = Un+i- this is because, as we remarked above, the boundary height around 
the boundary of Dn+i is approximately Un — {\/2)L?'^^. We deduce therefore that, with high 
probability, h^'^^{f) < Un+i = ipu^+iif) for t E [tn+i,L^] and a similar argument works for any 
other face / S {G')* . The claim at step n + 1 is then proven. 



9.3. Gaseous phases and entropic repulsion. Theorem 5.2 has been formulated for three 
specific - though quite natural - graphs. While, as explained in Remark 5.4, our method can be 
extended to a wider class of graphs, there is no hope it works for a general infinite bipartite graph. 
We would like to convince the reader there is a good reason for this. One of the main steps of our 



argument (cf. Section 8.3) is to prove that, given two height functions <h?^ after a step of the 



(fast) dynamics the mutual volume has not increased in average. The proof of Proposition 8.6 
shows that this is true (for the hexagonal, square and square-hexagon graphs), independently of 
the boundary conditions (in particular, for almost-planar conditions, independently of the slope 
(s,t)) and independently also of the presence of floor/ ceiling constraints. While obtaining the 
mixing time upper bound of order L^"*"*^ requires considerable extra work, the volume decrease 
result implies rather directly [13] that the mixing time is at most polynomial in L, since (i) the 
maximal volume difference between two configurations is of order and (ii) the ratio of mixing 
times of the fast and original dynamics is at most polynomial in L. 

Now take for instance the square-octagon graph, with almost-planar boundary conditions 



of slope (0,0). As we remarked in Section 3.2, in this situation the infinite- volume Gibbs 



measure /io,o corresponds to a "gaseous" (or "rigid") phase: the height h{f) at a face / has 
bounded variance and the random variables h{f),h{f') are essentially independent for /, /' far 
away. This is very reminiscent of the situation in the classical (2 + l)-dimensional Solid-on-Solid 
(SOS) interface model [21] at low temperature [Ij. Let us just recall that the SOS model 
describes an interface with integer heights 4'{x), labelled by a; G Z^, with measure proportional 
to Qwi~P'^\x-y\=i I'Pi^) ~ 4'{y)\) ^-iid zero boundary conditions (p = around a L x L box. 
Recently, it was proved in [3| that, when a floor at height zero is present (i.e. when heights are 
constrained to be non-negative) the mixing time of the Glauber dynamics for the SOS model is 
exponentially large in L. This effect is due to the rigidity of the interface and to the presence 
of the floor, which pushes the interface to a height of order logL (entropic repulsion [2]). If 
a "gaseous" phase does behave like a SOS interface, which is very likely, then that the dimer 
dynamics for the square-octagon graph with almost-planar b.c. of (0, 0) slope and floor at 
height zero has also a mixing time growing exponentially with L. As a consequence, for the 



reasons exposed above, the volume decrease cannot hold as stated in proposition 8.6 (i.e. for 
general floor/ceiling constraint and boundary condition) for the square-octagon graph (it cannot 
be true that the mutual volume between two arbitrary height functions decreases on average, 
under some reasonable "fast dynamics", for the square-octagon graph with general boundary 
conditions and floor /ceiling constraints). If on the other hand the volume decrease holds only 
for particular boundary conditions and floor/ceiling, then the mathematical mechanism for that 



must be considerably more subtle than in Proposition 8.6 



10. Mixing time lower bound 



In this section we establish the mixing time lower bound of Theorem 5.2 Thanks to Propo- 
8.3, it is enough to show that at time Tq = eL^ the asynchronous dynamics started from 



sition 



the maximal configuration is still at variation distance at least l/(2e) from equilibrium. 
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The strategy is the following. First, we define (for each of the three types of graph in question) 
a special (very non-planar) boundary condition p G il. and finite domain Wl C G of diameter 
of order L for which it is easy to prove that, starting from the maximal configuration, the drift 
of the volume is lower bounded by —cL. Therefore, after time To = eL^ the eroded volume 
is at most ceL^ and the configuration (call it Mtq) is still away from its (non-flat) equilibrium 
shape. Next, a monotonicity/coupling argument allows to deduce that, again at time Tq, the 
conflguration Mtq , evolving this time in our original domain G' with the almost-planar boundary 
condition m we are interested in, is above Mtq and that it is also far from its typical (flat, this 
time) equilibrium shape. 





Figure 8. "Pyramids" for the square, hexagon and square-hexagon graphs. The 
matchings are intended to extend over the whole inflnite graph. The dotted lines 
separate the domains tiled with different periodic matchings of extremal slope. 

10.1. Pyramids. For each of the three graphs (square, hexagon, square-hexagon) consider the 
special matching p £ ("p" for "pyramid" for reasons to become clear soon) of the whole G, 
defined through Fig. [8| Note that in all three cases G is divided into a finite number of infinite 
domains, separated by dotted lines (three domains for the hexagon graph and four for the square 
and square- hexagon graph): each domain corresponds to a vertex in the Newton polygon N of 
G and it is tiled with the periodic matching corresponding to that vertex (cf. Fig. [T]). The 
central face /o is assumed to contain the origin of and we fix the height there to 0. The 
large-scale height function 9 x i— t- H{x), obtained by rescaling the lattice spacing and the 
heights by 1/L — t- while keeping /o at the origin of the plane can be described as follows. For 
each vertex v of N, take the plane of the corresponding slope which contains the origin of M^, 
and let Sv be the half-space below it. The intersection of all the Sv with v ranging over the 
vertices of is clearly a pyramid 11 with vertex at the origin of M^. The boundary of 11 gives 
the height function H(-). It is possible to prove (but we will not need this directly) that the 
discrete height function associated to matching p is given by h{f) = —D{f, /o), cf. Section 4.1 
This observation could be used to build "pyramids" in a systematic way, for other graphs. 

Remark 10.1. For the hexagonal lattice, the pyramid p just corresponds (in terms of stepped 
surfaces, cf. Fig. to the surface of the corner of an infinite cube with vertex at the origin of 
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Figure 9. The subgraph Wl and the associated beads of configuration p. 
Threads are in Ught blue and white and beads are in sohd black. Arrows mark 
the orientation of threads. 



Next, we need to introduce a finite subgraph Wl C G, with the face /o in the center. This 
is defined through Fig. |9j for the hexagonal lattice this is the portion of G included in a 
(2L + 1) X (2L + 1) X (2L + 1) hexagon, for the square lattice it is a (2L + 1) x (2L + 1) 
square and for the square-hexagon lattice it is the portion of G contained in a suitable lozenge 
(delimited by a full line in the picture) whose diagonals contain 2L + 1 hexagons. Note that the 
maximal configuration in ^p,Wl is just p: indeed, just observe that none of the beads in Wl can 
be moved to a lower transversal edge in the same thread. 

Proposition 10.2. Consider the asynchronous fast dynamics on the finite graph Wl defined 
above, with "pyramid" boundary condition p. Let MJ^""^ (resp. M^'^ ) denote the state at 
time t of the dynamics started from the maximal configuration, which is nothing but p (resp. 
from the equilibrium uniform measure tTp^Wl) ^'^^ -^t be the sigma-algebra generated by 
{M^°-^ , Ms'^ , s < t}. Remark that M^'^ is stationary. Let Vt = X]/gWi,[^a/™^(/) — ^m^"'!/)]- 
Then, for every t ^ t' 

E[Vt - Vt,\Tt>] ^ -CL{t-t') 

for a certain constant G. 

Proof. The two evolutions can be coupled so that the one started from p always dominates the 
stationary one. Therefore, it is sufficient to prove that, given two height functions < h^, 



the initial-time drift of their volume difference, see (8.4), is lower bounded by —GL. Let d~WL 
denote the set of faces in Wl at graph-distance at most r from the exterior of Wl and assume 
the following 

Claim 10.3. There exists some finite r such that the initial-time volume drift is zero whenever 
h'^ and h^ differ at a single face f G Wl \ d~ Wl ■ 



As in the proof of Proposition 8.6, take a sequence of configurations /i(o)) • • • ' ^(fc) such that 
/i(o) = = h~ and is obtained by via a rotation that decreases the height at 

some face /. Along the sequence each face in d~WL can be rotated at most a number 

c(r) of times (since its height difference w.r.t. a face outside Wl at distance r can take at most 
c(r) = 2r + 1 values). Then write the volume difference between /i+ and h~ as a telescopic 
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sum of volume differences between h 
WlXO'Wl give zero drift (cf. Claim 
lower bounded by 



and h 



{i-iy 



the terms corresponding to a rotation in 



10.3) and those with rotation in d Wl give a contribution 



1 (cf. Remark 8.7). Then, one gets that the volume drift is lower bounded by 

□ 



-c(r) times the number of faces in d Wl, which is of order L, and Proposition 10.2 follows 



Proof of Claim 10.3. In order to follow, the reader should have in mind the proof of Proposition 
8.6: there we proved that the initial-time volume drift is zero except for "boundary effects", 



and here we show that indeed the boundary effects are not there sufficiently far away from 
the boundary. As in the proof of Proposition |8.6[ we consider for definiteness only the square- 
hexagon lattice, the other two cases being very similar. Note that there are 2L + 1 threads 
intersecting Wl (we label them 7^, i = —L, . . . ,L from left to right) and that thread 7^ contains 
L + 1 — |z| beads in Wl (we label them 6*-, j = 1, . . . , L + 1 — |i| from the lowest to the highest 
with respect to the orientation of the thread, see Fig. 10 a) for a schematic drawing). This 



geometric structure is essential for the proof of the Claim and, as the reader can check from Fig. 
[9J it is common to the "pyramids" of the three types of graphs. The effect of the boundary of 
Wl is simply that the beads of thread ji are constrained to stay strictly below some transversal 
edge ef E 7^ and strictly above some transversal edge e~ G 7^ (such edges are surrounded by 
circles in Fig. [Tol^a)). 

By symmetry, we assume that the face / belongs to 7^ for some —L + r<i<0 (the case 
—L < i < —L + r is excluded otherwise / G d~WL) and assume for definiteness that / is a 
hexagonal face (the argument is essentially identical when / is a square). Then, is obtained 
from h 

(Case 1) 



by moving a bead 6^one "transversal edge" lower along thread 7^ (with the notations 



of the proof of Proposition 



8.6 



b = 6*- moves from ei to 62). We need to show that, if 
/ G Wl \ d~WL, when threads of the opposite parity than i are updated while the threads with 
the same parity as i are frozen, thread 7j_i contributes exactly 1/2 to the change of volume (the 
same holds for 7i+i). 

We distinguish two cases, represented respectively in Figures [lO|b) and[To)(c): 

• l<j<L + l — Recall that from Section 

(which was called 



that, according to the position of 6* 



(which was called b' there), either 6* 



2 

-1 

has one more available 
transversal edge (called ey) when the configuration before the update is rather than 
h~ , or otherwise 6*~i (which was called b~) has one more available transversal edge 
(again ey) starting from h~ rather than from . Say that the former is the case. 
As discussed in Remark 8.7, the volume change due to thread 7i_i is 1/2 unless the 



boundary conditions prevent b~^ from moving down to 67, i.e. if 67 is not higher than 



this however cannot be the case, since 67 is clearly higher than b 



which is 



itself higher than e^_i 
j = L + 1 — \i\ (or j 



(here we use that j — 1 > 1, i.e. the bead is in Wl)- 
= 1, by symmetry). The argument is slightly different since this 



time neither b~^ = ^ nor b' = exist (since j>L + l — |i — 1|), or equivalently we 
can imagine that 6"^ , b' are higher than e 



1' 



ef_2 respectively (i.e. they are outside Wl)- 
Since the edge 63 of Fig. |6] is at distance of order 1 from /, we deduce that if / is at 
distance larger than some finite r from the boundary of Wl then b' is above 63. In this 
case, from Section 8.3, we get that 6^ has one more available position (transversal edge 
67) when the update of 7j_i is performed starting from configuration rather than . 
It is clear from Fig. |9] and [6] that, if / is at distance at least r, with r sufficiently large. 



from the boundary of Wl, the edge 67 is lower than 
b~ to actually reach it. 



-1' 



so there is no obstruction for 
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Figure 10. (a): A schematic view of the threads and beads in Wl (here, L = 3). 
(b): The case 1 < j < L+1 — |i|. Here, i = Oand j = 3. Boundary conditions only 
prevent beads from entering the gray region outside Wl, so they cannot prevent 
6"'" from moving down to ey. (c): The case j = L + 1 — If f £ d~WL then e-j 
can be at e^^, in which case it is not an ahowed bead position; otherwise, ey is 
in Wl (since it is at distance of order 1 from /) and b~ can reach such position. 



□ 

Volume 




Figure 11. An iUustration of the proof of the lower bound. The pyramid and, 
above it, the maximal configuration with planar boundary condition m. The top 
of the pyramid is 0{L) above the typical "almost-planar" equilibrium configura- 
tion (wiggled line). After a time eL^, the pyramid has only lost a volume eL^, so 
its top is still well above the "almost-planar" surface. By stochastic domination, 
the same holds for the maximal evolution with boundary condition m. 



10.2. Lower bound on mixing time. Recall Definition 2.2 the finite graph G' we are inter- 
ested in is the portion of G enclosed in Lf7, where C/ is a smooth bounded domain of M^, which 
without loss of generality we can assume to include the origin in its interior (so that /o is at 
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distance of order L from the boundary of G'). Now consider the finite graph Wg C G defined 



in Section 10.1, and choose £ to be the minimal integer such that Wi D G' (it is clear that 
£ = 0(1)). 

Consider first the evolution in Wg w ith boundary condition p (and started from the maximal 
configuration p): by Proposition 10.2 and Markov's inequality we have that at time To = eL^, 



nhTM<-e'/'L) = 0{e'/'') (10.1) 

(indeed observe that, if /i(/o) < —e^^'^L then the eroded volume is at least const x (e^/^L)^ ~ 
LT^er^^^ while the average eroded volume is of order LTq). 

Now recall that the boundary condition m we are interested in is almost-planar of slope (s, t): 
up to an irrelevant global change of the heights by an additive constant, we can assume that 
its height function hm{-) is within 0(1) from a plane of slope {s,t) and containing the point 
(0, 0, — 2e^/^L). If hp{-) is the height function associated to the "pyramid" matching p and if e 
is sufficiently small, hm{f) > hp{f) for every / along the boundary of G' (the height difference 
is actually of order L). To see this, just recall that the slope of the height function associated to 
p is negative and maximal (in absolute value) along any straight line starting from /q, while the 
slope (s, t) is in the interior of the Newton polygon, so it is not an extremal slope (see Figure 



11 ). Note that e has to be taken very small if the slope {s, t) is very close to dN. 

We have then that, by monotonicity, the evolution in G' , with boundary condition m and 
started from the maximal configuration, is stochastically higher than the restriction to G' of the 
evolution in We with boundary condition p. As a consequence, for the former dynamics we have 



that (10.1) still holds. Thanks to Theorem 3.7, at equilibrium (with almost-planar boundary 



condition m) the typical equilibrium height at /o is around — 2Le^/^ and the probability that 
it is higher than — Le^/^ tends to zero with L. This suffices to conclude that at time Tq the 
variation distance from equilibrium is still close to 1. □ 



Appendix A. Gaussian behavior of fluctuations 



In this section we prove Theorem 3.6 The proof is quite independent of the rest of the paper. 
The main tools are the results from [8j that were stated in Section [sj The proof is very similar 
to that of im Section 7] but the setting here is more general and we give a more explicit control 
of the "error terms". 

From the definition of height function, we know that h{f) — h{g) is determined by the dimers 
crossed by a path from / to g. In particular, for the square, hexagon and square-hexagon graphs, 
assume that / and g are on the same thread (cf. Section [7]): then, h{f) — h{g) is determined 
just by the number Nf^g of such dimers. Label ei, . . . , the transversal edges between / and 



g: it is easy to realize, using Theorem 3.2, that the set of occupied edges among the forms a 
determinantal point process: in other words, the probability ^s,t{^ii £ . . . , e^^ G M) can be 
written as det{A{ia,ib)i<a,b<ra) for a certain k x k matrix A directly related to K~l. Now, it 
turns out fTj that for the hexagonal lattice such matrix is Hermitian. In this case, a well-known 
theorem by Costin and Lebowitz (cf. for instance [19j) implies that is distributed like the 
sum of independent Bernoulli random variables, whose parameters are the eigenvalues of A. In 
particular, if the variance of Nf^g diverges as A; — >• oo, the variable [Nf^g — ^s,t{Nf^g)]/ ^yVar{Nf^g) 
tends to AA(0, 1). Unfortunately, for the square and square-hexagon graph the matrix A is not 
Hermitian and has complex eigenvalues, the Costin-Lebowitz theorem does not apply and the 
asymptotic moments have to be computed otherwise. 
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A.l. Choice of paths. Fix an integer k. We are interested in the behavior of 



M{f,f') 



(A.l) 



when the distance between / and /' goes to infinity. Remark that for any path C on G* from 
/ to /' that only crosses edges in the positive directiorj^ h{f') — h{f) is exactly the difference 
between the number of dimers crossed by the random matching minus those crossed by the 
reference matching, so M(/, /') = X]g[leeAf — Ms,t(6 £ M)] with the sum over all edges crossed 
by C. We will compute moments of higher order by taking k such paths Ci, . . . ,Ck ■ 



M{f,f') 



eiGCi 



(A.2) 



where we write 5e for leeM- 




Figure 12. The paths Ci, . . . ,Ck (here k = 5) along which is computed the 
height, displayed at large scale so that they look like continuous curves. They 
are macroscopically away from each other except near / and /' where they are 
periodic (so that on large scale they look linear). 



In principle we can choose the paths Ci . . . ,Ck freely as long as they only cross edges in the 
positive direction. In practice, for reasons that will be clear later we will adopt the following 
construction, illustrated in Fig. 12 Fix ri,r]' > 0. Inside balls of radius ijL around / and /', 



the Ci are portions of length rjL of infinite periodic paths (cf. Definition 2.1) and have mutually 
different asymptotic directions. Outside of these balls the paths stay at distance at least rj'L 
of each other and their length is of order L. Furthermore rj, r]' and the infinite periodic paths 
depend on k but not on L. 



A.2. Exact simplifications. Fix k edges ei = (wi,bi) G Ci, 
consider the corresponding term (written n(ei, . . . , e^)) of the sum (A.2). 



(wfc,bfc) G Cfe, and 



We write 6i for Se^, K^j for Kg /(bj,Wj) and Kij for i^s,f(wj, bj). We also write the set 



of permutations of {1, . . . , A;}, the set of permutations without fixed points and 5^ those of 
order 2. Given a S Sfc, we denote e{a) it signature. 

First we express each term n(ei, . . . , e^) of the sum (A.2) as a sum over permutations with 
no fixed points. 



^exercise: prove that such a path exists for every /, /' 
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Lemma A.l. Given edges ei, . . . , as above, we have 



71=1 n=l f^^Sk 

This is exactly like pi Lemma 21]. 

Now we replace K~^^^-^ by its asymptotic expression from Theorem 3.4 To lighten notations 

we use the following conventions. Recall that Cj = (wj, b^) and write (in a unique way) bj = bj + 
{xi,yi),Wi = Wj + {x[,yl) with bj, w, in the fundamental domain Gi and {xi,y.i) G 1?, {x'i->Vi) £ 
1? . Note that {xi — x[, yi — y[) is a vector of order 1. Then let Ui := Us,tii>i), Vj := Vs,t{'^j) and 
(f)i := (j){xi,yi), Xij =x[- Xj and yij = y'i - yj. We then get 



n(ei,...,e,)/(n^nn) = E ^(^)n 

1 



n=l 



n=l 



2'K{4>n - 4>a{n)) 



+ complex conjugate 



+0 



_|_ ^2 

n<7{n) ^na{n) 



E n(^".-(«) + <Hn) + Rn,ain))- (A.4) 



n=l 



We then expand the product over n and we call "dominant terms" those such that no R factors 
appear and, in addition, such that within each cycle of a either only A terms or only A* terms 
appear. All other terms of the expansion are called "error terms". We first show that in the 
dominant terms we can assume that all cycles of the permutation a are of order two. This comes 
from the following result: 

Lemma A. 2. If denotes the set of cyclic permutations of {1, . . . ,£}, then 

E n ^«-(") = E ^(^) n 



n=l 



n=l 



2^ 



•^ncr(n) yn(T{n) 



4>n - 4>a(n) 



na{ri) 



\ 



n=l 



if£>2. 

Proof. Recall that e(cr) is constant for cyclic permutations, thus it is enough to show that 



(A.5) 



En 



TT T r ■^nain) ynain) 

UnVa{n)Wo ' ' Zq 



- „ (Pn <Pa(n) 



as soon as £ > 2. The second equality is purely algebraic and is given in |1H Lemma 7.3]. □ 



Decompose a permutation a into cycles and remark that in equation (A.3) and (A.4), the 
sum over permutations with a given cycle structure can be factorized as a product over cycles. 
Then, Lemma A. 2 implies that the dominant terms in 

k 



n=l 
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exactly cancel each other so only "error terms" are left (recall that S"^ is the set of permutations 
without fixed points, and with only cycles of order 2). Altogether, we have proven (cf. ( A.4[ )) 

eiSCi efeeCfco-g52 \n=l / \n=l / 

In particular, if k is odd then there are only "error terms" because 5^ is empty, while, using 
equation (A.3) separately for all pairs, 

/i,,t(M(/, ff^) = ^ f^sAMif, f'ff + error terms (A.8) 

= 92k [Var^^ t {h{f) - h{f'))] ^ + error terms (A. 9) 

where g2k = {2k)\/{2^k\) is the Gaussian moment of order 2k. We will prove in the following 
section that the error terms are negligible, i.e. give a contribution to fj,s,t{M{f, f')^) that is 
much smaller than [Varf^^ ^{h{f) — h{f'))]^/'^. 

A.3. Controlling the "error terms". Recall that there are two kinds of "error terms" in the 
expression (A. 7) for the A;-th moment of the height fluctuation: those which contain both A and 
A* factors within the same cycle of a but but no R factor (recall ( A.4)), and those which contain 
R factors. The former will be shown to be "small" because they oscillate and give a negligible 
contribution when summed over ei, . . . , e^, while the latter are small because their denominator 
contains at least one factor — 4'(7{j) \ more than in the dominant terms. For simplicity purpose 
we assume in the following that a contains only a single cycle. The general case is essentially 
the same since an error term can be factorized as a product over the cycles of a. 
We will need the following elementary estimates: 

Lemma A.3. Define 



One has Fk{L) < LFk-i{L) and Fk{L) = 0((logL) 



\_k/2\' 



Proof. To get Fk{L) < LFk-i{L), writing Fk-i{L) either as a sum over di, . . . ,dk-i or over 
d2, ■ ■ ■ ,dk it is not difficult to see that 

2LFk-i{L) - 2Fk{L) = (di + ^2) • • • (4-2 + 4-i) (4 + 4-iK4 + + 4) - ° 

(use the fact that the denominator is symmetric under the exchange of di with dk)- As for 
F].{L) = 0((log L) L^/^-l ), this is an easy computation if A; = 2. By induction on /c, the proof is 
concluded if we show that Fk{2L)-Fk{L) = 0{kFk_2{L)) for /c > 4 and Fk{2L) - Fk{L) = 0(1) 
for A; = 3. To see this, one notes that the dominant contribution to Fk{2L) — Fk{L) comes 
from the terms where one of the variables is in [L, 2L\ while all the others are smaller than L: 
altogether, this gives 

^ 1 

, {di + d2)... {d^-i + {L + d,)){{L + di) + di+i) ...{dk + di) ^ ^^l?^k-i{L). 

For A; > 3 use Ffc_i(L) < LFk-2{L) and for A; = 3 note that F2{L) = 0{L). □ 



38 



BENOiT LASLIER AND FABIO LUCIO TONINELLI 



The first remark about the computation of iJ,s,t{M{f, f')^) is that we can restrict ourselves 
to the cases where either all edges are in the ball Bf{r]L) of radius "qL around either / or 
in the analogous ball around /'. Indeed, call di the distance of ei from / and observe that, 
since the map (j) is non-degenerate and the paths Ci are almost linear in Bf(7]L), one can bound 



— (pjl: from above and below, by a constant times {di + dj). Then the sum of (A. 4) with, for 
example, ei, . . . , ek-i in Bf{r]L) and out of Bf{r]L) is of order (1/ L)Fk-i{r]L) < Fk-2{'nL) = 
0((logL)L'=/2j-i) ^ [Far/,, - h{f'))]''l'^. Terms with more than one outside of Bf{r]L) 

are even smaller. To fix ideas, we will assume that all are in Bf{rjL). 

Consider now an "error term" containing a factor R (if it contains more than one, the argument 
is similar). Summing over e^, this gives a contribution of order 

^ , {di+d2)...{dk + di)di + d2 ^ {d2 + da) • • • {dk-i + dk) ^ cW ^^'"^^^ 

ai,...,(ife ^ r;L 0,2, ■■■,dk u=d2 

^ const X ^ , / '^''^f' =0{Fk-,{vL))<<[Var,^_Mf)-hU^^^^^ 

^ [d2 + (13) . . . [dk + 0.2} a2"/!c 

d2,...,dk 

We still have to deal with the error terms including no R factors but both A and A* within 
the same cycle of a. Recall that for simplicity of exposition that a is assumed to have a single 
cycle. Omitting the product of the factors ibi/(27r), these terms are of the form 

E e^['^(— )+^(-'--)]C(e„ . . . , e.) n 11 . \. ^ ^ (A.12) 

where ej runs over the edges crossed by path Cj and: 

• J (resp. J*^) is the set of indices n G {1, . . . , /c} for which we take A^a^n) (resp. 



in the expansion of the product (A.4). Note that both J and J'^ are non-empty, proper 
subsets of {1, ... , k}; 

• C(ei, . . . , Cfc) depends only on the types of the k edges (two edges being of the same type 
if they are related by a translation of Z^): 

• r, s are linear functions: if J' = {i G J : o"~^(i) G J*^} and J" = {i G J"^ : C7~^(i) G J} 
(remark that |J'| = \ J"\ 7^ otherwise o" would have more than one cycle) and writing 
Wo = exp{i6w), zo = exp{i6z), 

r{xi, ...,Xk) + 5(2/1, ...,yk) = 26l^(^ " X] ^'^^ + ^^^^X] ^'^ ~ X] 

aeJ' aeJ" aeJ' aeJ" 

Splitting the sum over the different types of edges, we can assume without loss of 

generality that all edges in each path are of the same type (so that C(ei, . . . ,efc) becomes a 
constant; edge types can be different in different paths) and that edges in the path j are obtained 

one from the other via translations by an integer multiple of some v^^^ = {v[^\v2^) G Z^. The 
edges in the j-th path will be labeled by an integer dj, which runs from 1 to Mj = 0{r}L) and 
r(xi, . . . , Xk) + s{yi, . . . , Uk) becomes is a linear function of di, . . . , dfc. 

Assume without loss of generality that 1 G J'. From the discussion above we obtain that 

(r + s)((ii,...,4) = {r + s){<d,d2, . . . ,dk) + 2{e^v^^H e,v^^^)di. 
Despite the fact that wq,zq are unit complex numbers with phase different from 0, vr it could 
happen that 6 := 2(6*^7;^^ + O^v^^^) is a multiple of 27r, so that exp(i(r -|- s)) is independent 



DYNAMICS OF PERFECT MATCHINGS 



39 



of di: this can however be avoided if the asymptotic directions and the period of path Ci in 
Bf{r]L) are chosen suitably (we skip tedious details on this point). 

We separate the sum over ei from the others and we make a summation by parts to get: 

gi(»-+«)Ul=o 

const X 2^ ^ = ^ (A. 13) 

62, 



'^2 - <P3) ■■■ {(t>k-l 



V-/ V- i0d\ A(^i((^fc + <^2 - 2(^i(di) - A(^i) / 1 



(Mdi) - (^2) {A - Mdi)) {Mdi + 1) - (t>2) {A - Mdi + 1)) \d2dk 



where the 0(...) comes from boundary terms in the summation by parts (note that e*®°' 
is bounded since / O(mod 2-k)), A0i = (pi^di + 1) — (/'i(di) is constant by linearity of cp and 
for simplicity of notation (pi can denote either (pi or its complex conjugate. We can thus bound 
\(pi — (pj\ by di + dj as before and (up to a constant factor) the absolute value of (A. 13) by 

I (y '^d, + d, + d, M^ (au) 

(d2 + d3)... (4-1 + 4) V4"(di + d2)2(di + 4)2 ^44^"/ ^ ^ 

_ 1 fsp (2^1 + 4 + 4)(4 + 4) ^ Q f d2 + dk 

~ (4 + 4).-.(4-i + 4)(4 + 4) (4 + 4)2(4+4)2 V c^24 

(A.15) 

= 0(Ffc„i(r?L)) « {Var^^^Xm " /^(/O)]'/' (A.16) 
because the parenthesis in the second line is clearly bounded. 
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